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COMPUTATIONAL TECHNIQUES FOR SOLAR 
WIND FLOWS PAST TERRESTRIAL PLANETS 
— THEORY AND COMPUTER PROGRAMS 

by Stephen S. Stahara, Denny S. Chaussee, 
Barbara C. Trudinger, and 
John R. Spreiter* 

Nielsen Engineering & Research, Inc. 


SUMMARY 


Theoretical analysis and the development of user-oriented computer 
programs were carried out for the purpose of developing computational 
techniques for predicting the interaction of the solar wind with terres- 
trial planets. The procedures are based on a single-fluid, steady, 
dissipationless, magnetohydrodynamic model and are appropriate for the 
calculation of axisymmetric , supersonic, super-Alfvenic solar wind flow 
past both magnetic and nonmagnetic planets. The actual calculations are 
implemented by an assemblage of computer codes organized into one program. 
These include finite-difference codes which determine the gasdynamic 
solution, together with a variety of special-purpose output codes for 
determining and automatically plotting both flow field and magnetic field 
results. , 

Theoretical results obtained using these procedures are reported for 
a variety of solar wind conditions and different obstacle shapes. 
Comparisons are made with previous results, and new results are presented 
for a number of solar wind flows. The computational programs developed 
under this work have been documented and are presented in a general user's 
manual included as part of this report. 

INTRODUCTION 


The gasdynamic and magnetohydrodynamic calculations of solar wind 
flow around magnetic and nonmagnetic planets that Spreiter et al . (refs. 
1-8) carried out some years ago have been and continue to be widely used 
in the interpretation of data measured in space around the Earth and other 
planets. The objective at the time those calculations were made, was to 


*Professor, Department of Applied Mechanics, Stanford University. 
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provide some theoretical results based on a fluid rather than particle 
description of the flow that might be compared with measurements in space 
to determine the applicability of such a description. During the inter- 
vening years, the usefulness and accuracy of those results have become so 
well established that they are currently being used, for instance, as one 
of the ways in which the magnetic field of a planet such as Mercury or 
Mars is estimated from measurements of fly-by or orbiting spacecraft. 

In such applications , the previous calculations leave much to be 
desired. They have been carried out for only a limited set of conditions 
such as obstacle shape, oncoming Mach number, angle between free-stream 
flow and magnetic vectors, etc., and are presented in archival publica- 
tions only in the form of plots from which results for other conditions 
must be determined by scaling and interpolation. The usefulness of the 
present theoretical model, therefore, would be considerably enhanced if 
similar but more detailed results for the specific conditions measured by 
a spacecraft would be provided by the application of a documented, user- 
oriented, economical computer code. 

That objective has been successfully accomplished in the work 
reported here. Based upon the identical theoretical model employed 
previously (refs. 1-8) , and discussed in more detail below, a new more 
versatile means of solution far superior in efficiency and generality to 
the former method is now available. Current advances in computational 
methods have been incorporated in these programs and, in fact, an entirely 
new blunt body code based on a recently published algorithm was developed 
to determine the flow field in the region surrounding the nose of the 
obstacle. The new methods confirm previous results and are readily 
applicable to additional cases involving different geometries and flow 
conditions. The entire procedure is fully automated and provides detailed 
flow field and magnetic field properties in a convenient output format. 

An automatic plotting capability for generating report-quality plots is 
also included. 
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internal energy, eq. (5) 
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entropy 

incremental distance along streamline 
time 

velocity components associated with the (X,R,u>) coordinate 
directions, respectively 

column matrix defined by eq. (29) 

column matrix associated with IMP code, equal to U/J 
column matrix defined by eq. (34) 
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P 

R 

t 

0 
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parallel, perpendicular, and normal magnetic field components 
as defined in eq. (42) 







ANALYSIS 


General Considerations 

An overall description of solar wind flow past terrestrial planets, 
including an account of the principal features of the interplanetary solar 
wind and a description of the physical basis of the continuum fluid model 
are provided in references 3 and 6. Although details of those accounts 
won't be repeated here, several specific points relevant to the present 
work will be reviewed. 

The fundamental assumption underlying the theoretical analysis of 
large-scale features of the interaction of the solar wind with terrestrial 
planets is that the average bulk properties of the flow can be described 
adequately by the continuum equations of magnetohydrodynamics for a 
single-component, perfect, dissipationless (perfect electrically con- 
ducting, inviscid, nonheat-conducting) gas. Theoretical justification of 
this point has not yet been established, and proof remains essentially 
qualitative at present. The primary justification for use of the 
continuum fluid model is the outstanding agreement of the qualitative 
results predicted on this basis with those actually measured in space. It 
appears that the continuum model is capable of accounting both for many of 
the details as well as the broad features of the observations. 

The primary emphasis on development of the continuum model has been 
on steady flow past the Earth (refs. 3, 4) , a planet with a strong 
magnetic field. Initial applications to nonmagnetic planets have also 
been made (ref. 5). These particular applications are the ones toward 
which this study is directed. 

The following subsection provides the mathematical basis of the 
interaction of the solar wind with planetary magneto/ionospheres. The 
governing differential equations as well as relations for discontinuities 
present in the flow are discussed. The next two subsections provide the 
determination of the magnetosphere and ionosphere boundary shapes. The 
following two subsections discuss the calculation of the gasdynamic and 
magnetic field properties. 
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Mathematical Representation of Solar Wind- 
Magneto/Ionosphere Interaction 


The differential equations governing the motion of the solar wind 
are the magnetohydrodynamic equations for steady flow of a dissipationless 
perfect gas and can be written as follows 


y.pv = 0 

V 'Vi 



p (V*V) V + Vp 

'Vi V 'Vi v 

Curl (B x v) = 

<\i 'Vi *\i 

div B = 0 

v 

(V*V)S = 0 

'V/ ^ 

S - S = C in 
o v 


-(1/4tt) B 
% 

0 


P/Pp 

(p/p Q ) Y 


X 


Curl B 

v 'Vi 


J 


and apply in the region exterior to the magneto/ionosphere boundary, 
shown in the sketch below. 


(1) 


as 


Solar 

Wind 



In these equations , p, p, S, and V refer to the density, pressure, 

entropy, and velocity of the gas, B denotes the magnetic field, and 

% 

y = Cp/C v is the ratio of specific heats, equal to 5/3 for a monatomic 
gas. The temperature T is related to the pressure and density by the 
equation of state for a perfect gas 


p = p R T/p 


( 2 ) 
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in which R = (C - C )y = 8.315xl0 7 ergs/gm°K is the universal gas 

p V 

constant, and y is the mean molecular mass nondimensionalized so that 
y = 1.6 for atomic oxygen. For fully ionized hydrogen,, y is .thus 1/2. 

The magnetic field B is expressed in terms of gaussian units. 

^ ‘ . • .. .* * 

We note that it is fully equivalent and convenient for some purposes 
to replace the entropy equation (V*V)S = 0 in equation (1) by the fol- 
lowing energy equation: 


div[pV(2 V 2 


+ h) + 


VB 

a. 


! /4tt - 


(B • V) B/4tt ] = 

'X, O/ % ■ 


(3) 


Because of the omission of dissipative terms in equation (1) , 
surfaces of discontinuity, such as the bow shock and the magneto/iono- 
sphere boundary indicated in the sketch above, may develop in the flow. 
Across these surfaces, continuous solutions of the dissipationless 
differential equations cease to exist. The flow is no longer governed 


solely by the differential equations (1) , but must . be supplemented by 
additional considerations. Mass, momentum, magnetic flux , -and^energy, . 
must be conserved across these surfaces, and these conditions lead to the 
following relations which relate quantities on the two sides of any such 
discontinuity: 


[pV | - 0 / „ 

n ; - , ... 

. [pV- n V; -+’ (p + .B 2 /8 tt ) n - B n B t v(4n] = 0 


[B nXt - StV = 
[B n> * 0 



[pV (h + V 2 /2) + V B 2 /4 tt - (B*V)B /Aw] = 0 I 
n n a. % n 


(4) 


Here, (n, t) denote unit vectors normal and tangential to the discontinuity 
surface, as sketched below. 
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The basis for important simplifying approximations to the governing 
tangential discontinuity conditions 


V = B = [p + B 2 /87 t] = 0 
n n 

[V. ] + 0; [B ] * 0; [p] + 0 


(9) 


which apply at the Earth's magnetosphere boundary are that the gas 
pressure p is much less than the magnetic pressure B 2 /8tt in the confined 
region interior to the magnetopause, and that p is much greater than 
B 2 /8tt in the exterior flow region around the forward part of the magneto- 
pause. Consequently, the boundary condition to be satisfied by the 
exterior flow at the magnetopause can be satisfactorily approximated by 
the limiting case of the tangential discontinuity conditions in equation 
(9) in which there is a vacuum 1 (p = p = 0) on one side, and no magnetic 
field on the other. Thus, the condition [p + B 2 /8tt] = 0 becomes 


■- <p>«0« no) 

while the equations governing the magnetic field interior to the magneto- 
pause are ' 

• » f • . 

div B = 0, Curl B = 0 (11) 

^ v % 


Furthermore 
the flow on the 
in equation (10) 


, it has been shown (refs. 3, 6) that the gas pressure of 
forward portion , of the magnetosphere boundary • — Op) 

— can be approximated . adequately by the Newtonian formula 


P = P st cos 2 ip (12) 


Here, ip is the angle between the outward normal to the magnetosphere 
boundary and the flow direction- of the. oncoming undisturbed solar wind, 
and p. is the stagnation or ram pressure exerted on the nose of the 
magnetosphere and is given by 

P S £ = K ' - ■ .(13) 

In this relation, K is a constant usually taken as one, but whose actual 
value is 
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For the high Mach number flows typical of solar wind conditions, K 
approaches 0.844 for y = 2 and 0.881 for y = 5/3. The important implica- 
tion associated with the introduction of the Newtonian approximation is 
that the calculation of the shape of the magnetosphere boundary decouples 
from the calculation of the external flow. 

•Implementation of the above analysis' requires a representation of the 
Earth's magnetic field; This is provided with sufficient accuracy for the 
present application by a magnetic dipole located at the center of the 
Earth and given by • 

B = -B (a /r) 3 (sin 8 e a + 2 cos 0 e ) (15) 

v eq e 0 r 

where B = 0.312 gauss is the strength of the field at the geomagnetic 
equator, a g = 6. 37x10 8 cm is the radius of the Earth, r is the geocentric 
distance, and 0 is the polar angle' measured with respect to the north 
geomagnetic pole which is at 78.6° north latitude and 70.1° west longitude. 
The problem for determining the magnetopause and the distortion of the 
confined geomagnetic field is then- reduced to finding a solution of 
equations (11) which has the dipole singularity equation (15) at the 
origin, and which satisfies the tangential- discontinuity conditions 
B n =. 0 and - B 2 /8tt" = Kp^V^ 2 cos 2 ip at- the unknown location of the magneto- 
pause. This problem is identical to the classical steady-state Chapman- 
Ferraro problem (ref. 10) specified over four decades ago, and based on 
particle rather than continuum concepts. Even at this stage, however, the 
problem still remains sufficiently complex that only a few numerical 
solutions have been carried out (refs. 11—13) at the present time. An 
early approximation which has proven quite successful in providing 
accurate magnetopause shapes (refs.' 9, 12, 14) is to replace the B n = 0 

condition with the condition that the magnetic-field intensity at the 
magnetosphere boundary is equal to 2f times the tangential component of 
the Earth's dipole field, where f is a constant usually taken as unity 
(refs. 15-17) . 
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The final mathematical statement for the free-boundary problem for 
the. shape of the magnetosphere boundary, then,, is summarized in the sketch 
below : 



This formulation leads to the following partial differential equation 
for the geocentric distance r to the magnetopause for the case in which 
the dipole axis is perpendicular to the solar-wind flow (ref. 18) 


1+3 cos 2 0 

1 

3R 

m 

2 


sin 0 

, 2 cos 0 

dR 

m 

2 

R b 

m 

R sin 0 
m 

9<j> 

> 


+ R^. 
m 

R 

3<j> 

. J 



sin 4> sin 0 - 


sm 


R 


D SR 

cos 8 m 
36 


cos 


3R 


m 


m 


R sin 0 3d) 

m 


(16) 


, , 7 1 / 6 

where R = r/D, D = a (f B /2 ttKp V ) is the geocentric distance to 
m e eq °° 00 

the magnetosphere nose, <J> is the geomagnetic longitude measured with 
respect to a line through the origin that is normal to both the dipole 
axis and the free-stream direction, and is equal to tt/ 2 when directed in 
the upstream direction, and 3tt/2 when directed downstream. 


A considerable simplification occurs if attention is confined to 
either the geomagnetic equatorial plane or the principal meridian plane 
defined by the dipole axis and the direction of the incident solar wind. 
These equations are; in the equatorial plane (0 = tt/2, SR^/30 = 0) 


14 



dR ™ 

^rr = R 

d<p m 


R 6 sin d) cos tb + /R® - 1 
m . m ■ 


cos 2 it - r 

m T 


m 


Tl/2 < 4> < 3 7T 
2 


8 R 

and in the principal meridian plane (<S> = ± tt/2, — = 0) 


dR 


R 3 + 1 


df = (R m tan 6) FTI 

m.. 


- tt/2 < 6- £ j 


( 17 ) 


(18) 


A sketch of the traces of the Magnetosphere boundary in the equator- 
ial and principal meridian planes is provided below together with the 
spherical (r, 9, <j>) coordinate system. We note in particular that the 




shape of the equatorial trace is a smooth curve, while that of the 
principal meridian contains a pronounced dent at the location of the 
neutral point where the magnetic field vanishes in the ideal theory. In 
all of the reported work to date (refs. 1-4, 6, 8) concerning flow fields 
past the magnetosphere, the magnetopause has been approximated by an 
axisymmetric shape obtained by rotating the equatorial trace of the 
boundary about the longitudinal axis. That this provides a very good 
approximation to the actual three-dimensional shape is shown in figure 1 
which displays in more detail the comparison of the equatorial and 
principal meridian traces of the magnetopause as calculated by equations 
(17) and (18) , respectively. Also shown on the principal meridian trace 
is a dashed line which forms a tangent surface across the cusped neutral 
point region, and which represents the free-streamline surface which will 
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most likely appear in nature (ref: 19).. We note that all three of these 
shapes are quite close to one another so that the rotated equatorial trace 
provides a satisfactory approximation to the three-dimensional shape. 
Furthermore, it has been shown (ref. 3, 13) that higher-order approximations 
to the Chapman-Ferraro problem for the magnetosphere shape are in reason- 
able agreement with the approximate shape given by the solution to 
equation (16) . 

Consequently, the boundary shape of the magnetosphere for which the 

majority of results are presented here is the rotated equatorial trace 

given by equation (17) . That shape has been incorporated in the associated 

computer codes developed herein as the default shape for magnetosphere 

calculations . .Tabulated ordinates (Y /D - vs. - x/D) for the equatorial 

m 

trace are provided .in . Table • 1 where is the cylindrical radial coordinate 
of the profile. - ■. ‘ i ; • ■ \ 

Nonmagnetic Planet —Determination 
of the Ionosphere Boundary 

. ' \ e' * 

The determination of the ionosphere boundary initiates from the 
assumptions that the ionosphere, or at least the outer part of it that 
participates in the interaction with the solar wind, is idealized as a 
spherically symmetric and hydrostatically supported plasma having infinite 
electrical conductivity, effectively bound to the planet and incapable of 
mixing with the solar wind, as indicated in the sketch below: 



16 ' 



This interior plasma is separated from; the flowing solar plasma by a 
tangential discontinuity across which the same, relations, given previously 
by equations (9) for the analogous problem for a magnetic ■ planet , must 
hold; . The assumption of hydrostatic support implies a quiescent iono- 
sphere, in which the bulk motions of-. the gas with ' respect to the planet are 
sufficiently small (V = 0) that equilibrium exists between the pressure 
gradient and gravity, viz. 


dp _ _ 


=? ■ -pg 


where p and p are the gas- pressure and density, r is the radial distance 

measured from the center, of the planet, and g is the acceleration due to 

gravity.. The variation of g is inversely proportional to r 2 ',- so that 

g = g (r /r) 2 where the subscript s denotes values at the surface of the 
s s 

planet. Since the density p is related to the pressure according to the 
perfect gas law equation (2) ,, equation. (19) can be. integrated to yield 


r dr 
exp -/ — 


where p is the pressure at some reference radius R n and H is . the , local 
scale height of the atmosphere given by H = RT/ g . If H is regarded as 
constant; that is, if variations of g and T with r- are neglected , equation 

■ • • * J ’ t • 

(20) can be integrated directly to yield 


P = P R .exp|- 


r - R. 


In view of uncertainties associated with . measurements of the atmospheric 
properties of Venus and Mars, the variation of p with r as given by equa- 
tion (21) was adopted in previous solar wind/ionosphere applications 
(ref. 5) and has also been. used herein. 

In the present application, it is also evident that the gas pressure 
p is much larger than the magnetic pressure B' 2 /8 tt. on both sides of the 
ionopause. Therefore, the discontinuity pressure balance relation 
[p + B 2 /8 tt] = 0 of equation (9) reduces to a simple equality between the 
ionospheric pressure given by equation (21) and the static pressure of the 
flowing solar plasma adjacent to the ionopause. Introducing as in the 
magnetosphere application the Newtonian approximation equation (12) for 
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the static pressure of the flowing plasma on the exterior boundary of the 
ionosphere, we arrive at the following equation for the pressure balance 
at the ionopause : 


p . cos' 
r st 


i|) = Kp V ' 

r oo oo 


cos 


^ = P R exp 


R . - R_ 
l R 


( 22 ) 


It is convenient to choose as the reference radius and location the stag- 
nation point on the ionopause; that is, R R = R q where R^ is the distance; 
from the center of the planet to the nose of the ionopause. This implies 
that 


p„ = p = Kp V ' 

and that at all points along the ionosphere boundary 


cos 2 ip = exp 


R. - R 
i o 

H 


( 23 ) 


( 24 ) 


The final mathematical statement of the free-boundary problem for 
c.Btermining the shape of the ionosphere boundary then is summarized in the 
s etch below: 



By relating the local angle ip to the slope of the ionopause, the following 
ordinary differential equation results for the ordinates of the ionosphere 
boundary 
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(25) 


dR. 
l 

d 3 


sin 2 3 - 2 /A - K 2 
2 (A - sin 2 3l 


where 

A = exp (- ( R ± - R q )/H] (26) 

and 3 is the angle measured from the subsolar point as indicated above. 
Results for various ionopause shapes obtained by integrating equation (25) 
for different values H/R q in the range 0.001 <_ H/R q <_' 1 . 0 are provided in 
figure 2. We note that the range of interest for planetary applications 
to Venus and Mars appears to be 0.01 < H/R < 0.30 (ref. 5). Tabulated 
ordinates (Y^/R^ ~ vs. - X/R q are provided in Table 1 for H/R q = 0.01, 

0.1, 0.2, and 0.5, where ib = R^ sin 3 is the cylindrical radial coordinate 
of the ionopause profile. 

Calculation of the Gasdynamic 
Flow Properties 

Computation of the gasdynamic flow properties consists of determining 
solutions to the differential equations and discontinuity conservation 
equations given by equations (7) . Because both the downstream tail region 
(far field) as well as the region in the vicinity of the obstacle nose 
(near field) are generally of interest in solar wind-magneto/ionopause 
interactions, the computational methods selected must be capable of 
determining this entire flow field. In view of the need to carry the flow 
calculation to an arbitrary downstream distance, the most computationally 
efficient procedure to employ is to subdivide the flow field into two 
regions, as indicated in the sketch below: 


Blunt Marching 
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Illustrated here are the essential features of the high supersonic Mach 
number flow typical of solar wind flows past terrestrial planets. Of 
special note is the embedded subsonic pocket, located at the nose of the 
magneto/ionopause. The presence of this pocket necessitates use of a 
computational technique capable of treating a mixed subsonic/supersonic 
flow. Downstream of this region, the flow becomes supersonic and remains 
so for the convex shapes typical of magneto/ionopause boundaries. In that 
ir field region, a more economical procedure than that used near the nose 
can be employed. 

We note that this subdivision of flow field and use of different 
solution procedures in each region is common practice for calculating 
hypersonic blunt-body flows, and was employed in the former solar wind 
method (ref. 3) as well as in a more recent application to space shuttle 
re-entry flows (ref. 20) . Indicated in figure 3 is a comparison of the 
former and present computational procedures which illustrates the break- . 
down of the various solution domains. The previous method treated the 
nose region by the inverse iteration method, and the remainder of the flow 
field by the method of characteristics. In light of recent advances, both 
of those techniques, particularly the inverse method, are now considered 
dated and inferior to more current methods. In the present study, those 
two methods have been superseded by: (1) a new axisymmetric implicit 

unsteady Euler equation solver (IMP) specifically developed under this 
contract (ref. 21) and which determines the steady state solution in the 
nose region by a time-marching procedure, and (2) a shock -capturing 
marching solution (SCT) which spatially advances the solution downstream 
as far as required by solving the steady Euler equations. As indicated in 
figure 3, the implicit unsteady Euler equation (IMP) code is used to 
determine the flow in the region between the stagnation streamline on the 
symmetry axis and a body-axis normal plane (x = constant) conveniently 
chosen at the 3 = 90° (x = 0) location. For the magneto'/ionopause shapes 
considered here, purely supersonic flow exists on this plane and the 
conditions there provide the starting data for the SCT code. The SCT code 
then marches the solution downstream in axis normal planes to the final 
specified location. 

Obstacle nose solution — implicit unsteady Euler equation (IMP) code .- 
The partial differential equations employed in the implicit (IMP) code are 
the unsteady axisymmetric versions of the gasdynamic Euler equations given 
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by equations (7). To implement the' calculation, those equations are 
written in conservation-law form under the generalized independent variable 
transformation f ' 


T * = 'T • ■' 1 - ' 

£ = £(T,X,R) • ' 

.n =■ n.(T f x,R). ^ . 


(27) 


as follows: 
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(29) 


J . ? X n R. ” .^R n X •. (30) 

In equations (27) to (30) , T denotes time, X is the axial downstream 
direction, and R is the cylindrical radia 1' distance ; p represents the 
pressure; p the density; u and v' the' velocity components in the X and R 
directions, respectively; and e the total energy per unit volume. The 
following equation relates the pressure, density, and velocity components 
to the energy for an ideal gas 

e = p/(y - 1) + p(u. 2 + v 2 )/2 (31) 
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Using the independent variable transformation equation (27) , the 
physical plane with boundaries 1-4 in X,Y,T space is mapped into the 
rectangular computational plane fj , n, and t space as shown in the sketch 
below. Generally, at each time step the transformation metrics are. not 
known beforehand and must be determined numerically. Integration step 
size is established by using the eigenvalues of the Jacobian matrices A 
and B, where A = 3E/3U, B = 3F/3U and U = U/J , E = U T U + 5 X E + £ R F)/J, - 
and F = (n T u . + H X E + n R F)/J. 



Physical plane Computational plane 


The boundaries, as indicated in the sketch, at which appropriate 
conditions are imposed on the present problem consist of: (1) the bow 

shock surface at which the Rankine-Hugoniot relations are satisfied, 

(2) the downstream outflow boundary where the flow is assumed to be 
entirely supersonic, (3) the obstacle surface at which a no flow condition 
in the normal direction is imposed, and (4) the stagnation streamline 
where symmetry conditions are implemented. Initial flow field conditions 
are determined by guessing a bow shock shape for the particular blunt 
obstacle at hand and by prescribing a Newtonian pressure distribution on 
the body. Since the maximum entropy streamline wets the body , that fact 
plus the known surface flow direction serve to determine the remainder of 
the flow properties on the obstacle surface. A linear variation for the 
flow properties between the bow shock and the obstacle is then prescribed. 
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This provides the initial flow field which is then integrated in a time- 
asymptotic fashion until the steady-state solution is obtained. 

The basic numerical algorithm used in the IMP code was developed by 
Beam and Warming (ref. 22) and is second-order accurate, noniterative, and 
spatially factored. In particular, the "delta form" with Euler time 
differencing is employed. When applied to equation (28) the algorithm 
assumes the form 

(I + Ax6 A* 1 ) (I + ATS n B n ) (U n+1 - U n ) = -Ax(S_E n + 6^F n + G) (32) 

where A and B are the Jacobian matrices, I is the identity matrix, 6 and 

n+ i ^ 

6^ are second-order', central difference operators, U = U(nAx) and At is 
the integration step size. 

Equation (32) is solved at the interior points only. It requires two 
4x4 block tridiagonal inversions at each time step of the integration. 

The solution procedure is as follows: 

n+ 1 

. 1. Define AU = U - U . 

2. Form the right-hand side of equation (32.) and store results in 

, , r,n+ i 

the U array. 

3. Define U = (I + Ax6 B ) AU and solve the matrix equation 

— n — ~ n + 1 — ^ ^n+ i 

(I + At6j.A ) U = U for U storing the result in the U array. 

— FI AnX 1 

4. Solve the matrix equation (I + Ax6' B ) AU = U for AU. 

n 

. ~n + i i ^ 

5. Obtain the values of U from the relation U = AU + U . 

6. Obtain the final solution after applying the smoothing; U n+1 = 
U n+1 - (c/8) S/J. 

A fourth-order smoothing ■ term S is usually required to eliminate 
nonlinear instabilities which may arise since the use of central differ- 
ences in the spatial directions results in a neutrally stable algorithm. 

Downstream solution — shock capturing steady Euler equations (SCT) 
code . - Since the shock-capturing technique employed has been established 
previously (refs. 23, 24, 25) only a brief description will be provided 
here. The steady, full three-dimensional Euler equations , are written in 
cylindrical coordinates in the form: 
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Here p and p represent dimensional pressure and density and (u,v,w) denote 
velocity components in the cylindrical coordinate directions (X,R,io) , 
respectively, where X is the axial downstream direction, R is the cylin- 
drical radius in an axis-normal plane, and to is the azimuthal angle 
measured in the axis-normal plane. This set of equations is closed by the 
use of energy conservation as provided by the equation for total enthalpy 


= h(p,p) + SL- = constant 


(35) 


where 

q = /u 2 + v 2 + w 2 (36) 

is the magnitude of the velocity vector and h(p,p) is the caloric equation 
of state . • 

The equations are then transformed to a computational plane by using 
the independent variable transformation 


X = X 


(R - R ) 

5(X,R,w) = (37) 

£ (to) = f (w) 
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where the subscripts (b,s) denote body and shock surfaces, respectively. 
This change of variable is employed in equation (33) and the result then 
rearranged into the following conservation law form (refs. 23,24,25). 


3U 3F 3G 

3X 3£ 3n 


H 


0 


(38) 


The finite-difference formulation of MacCormack (ref. 26) is used to 
integrate equation (38) with respect to the hyperbolic coordinate X to 
yield values of the conservative variable U. This procedure is calculable 
as long as the flow velocity in the axial (X) direction is greater than 
the local speed of sound. We note that for the axisymmetric obstacle 
shapes considered here, 3/3w = 0; however, the version of the SCT code 
employed has the capability for treating full three-dimensional geometries. 


Streamline and contour calculations .- A number of special purpose 
subroutines were written to calculate streamlines and contour maps of 
various flow field properties, and also to provide automated plots of 
these quantities. The streamlines are determined by integrating trajec- 
tories through the known velocity field as this procedure was found to be 
more accurate than a mass flow calculation. The calculation of a particu- 
lar streamline is initiated at the point where the streamline crosses the 
bow shock, as shown below: 



At that point, exact values of the streamline slope dR g /dX are known in 
terms of the local shock angle 6 g and free-stream quantities according to 
the relation 
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which is implicit in the blunt-body (IMP) and marching (SCT) code solu- 
tions. At other points in the flow field, the local streamline slope is 
given by the ratio of radial to downstream velocity, i.e. 



(40) 


and the streamline determination is made by stepwise integration in X of 
equation (40) according to a modified third-order Euler predictor/corrector 
method. Two-dimensional linear interpolation from the flow field grid 
points is employed to obtain the velocity components (u,v) required at the 
stepwise points along the streamline trajectory. Separate streamline 
calculations are made for the near field (IMP results) and far field (SCT 
results) because of the different coordinate systems employed in those two 
regions . 

Contours of various flow field quantities are determined for the 
entire flow field (near field plus far field) by using a modified version 
of the contour package developed by R. Sorensen at NASA/Ames Research 
Center. Pertinent details of that package are provided in the user's 
manual section of this report. 

Calculation of the Magnetic Field 


With the flow properties known from the gasdynamic calculations, 

determination of the magnetic field B proceeds by integrating equations 

a. ' 

(8) . Those equations are commonly interpreted as indicating the field 
lines move with the fluid. This leads to a straightforward, although 
tedious, calculation in which the vector distance from each point on an 
arbitrarily selected fluid line to its corresponding point on an adjacent 
field line in the downstream direction is determined by numerically 
integrating / v dt over a fixed time interval 6t. While that calculation 
determines the location of the field lines, the intensity of the magnetic 
field at any point may then be calculated from the relation 
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( 41 ) 
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where A H is the vector length of a small element of a flux tube. The 
sketch below clarifies these quantities for the case when the oncoming B 
field is perpendicular to the flow: 



where the open symbol (2) denotes locations of points on the streamlines 
corresponding to a fixed time interval At = AS^/V^. 

Because of the axisymmetry of the gasdynamic quantities and the 

linearity of the equations for B, Alksne and Webster (ref. 27) have shown 

a. 

that the magnetic field at any point in the flow field can be calculated 
by vectorially summing the contributions of the three component fields 
indicated in figure 4. At any point P, the magnetic field B is given by 
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where the subscripts ( " , , n) refer respectively to the contributions 

associated with the components of B m parallel to perpendicular to 

in the plane that contains the point P, the center of the planet, and the 
vector V^; and normal to the latter plane. 

For the contribution associated with the parallel component B , it 

00 1 I 

has been shown (refs. 28,29) that (Bp) it is proportional to pV throughout 
the entire flow field. Consequently, this component can be calculated 
directly from the gasdynamic solution by the expression 
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Similarly, the normal component (B /B ) can be shown to be equal to 
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( 44 ) 


where Rp is the radial cylindrical coordinate of the streamline, as indi- 
cated in figure 4. 


The determination of the contribution (Bp/B^)^ cannot be made directly 

from the gasdynamic properties as can the other two components. For this 

component, recourse must be made to the general computation associated 

with equation (41) . Details of the calculation procedure for (B,/B ) are 

<ViP 00 

provided in the user's manual section of this report. 
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RESULTS AND DISCUSSION 


In order both to verify the correctness of the procedures developed 
under this study as well as to demonstrate their flexibility and power for 
calculating solar wind flows for a variety of conditions, a large number 
of test cases have been run. A representative sample of the results 
obtained from those calculations is reported here. Because the implicit 
blunt-body procedure (IMP code) is new, that code was tested and compared 
to previously established theoretical methods to substantiate its validity 
prior to solar wind applications. We have prepared figure 5 to exhibit 
such a comparison. Here, flow field results predicted by the present 
implicit method are compared with 'those of other theoretical techniques 
and also experimental data for supersonic flow past a sphere with M^ = 

4.926 and y = 1.4. The variation of density and Mach number along the 
stagnation streamline are provided in the two plots on the left, while the 
variation of surface pressure and shock standoff distance with spherical 
angle are given in the two plots on the right. In all of these compari- 
sons, the present theoretical results are in excellent agreement with both 
established theoretical methods and experiment. 

Figure 6 exhibits a comparison of results predicted by the implicit 
method for a solar wind calculation. The upper plot of the figure displays 
the locations of the bow wave and sonic line for flow past the equatorial 
trace of the magnetopause for M^ = 8 and y = 5/3. The density distribution 
along the magnetosphere boundary and along the shock wave for this flow 
are given in the lower plot. Included in both of these figures is the 
solution calculated by Spreiter et al. (ref. 3) using the inverse method 
of Inouye and Lomax (ref. 30). Excellent agreement is obtained between 
the present implicit technique and the inverse method. Figure 7 provides 
analogous results for the equatorial trace of the magnetopause for a 
variety of Mach numbers typical of solar wind flows, and illustrates the 
Mach number capability of the present IMP code. 

With regard to solar wind flows past nonmagnetic planets, figure 8 
displays results for the bow shock location for M^ = 8 and y = 5/3 flow 
past various ionopause shapes with H/R^ = 0.01, 0.1, 0.2, 0.25, 0.50, 0.75 
and 1.0. Not included in that figure for clarity of presentation are the 
results of Spreiter et al . (ref. 5) obtained by using the inverse method. 
Those results and the predictions of the present implicit method given in 
figure 8 are essentially indistinguishable. 
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As an illustration of the geometric flexibility of the present 
implicit method and also as a critical test of its ability to capture 
embedded shock waves besides the bow shock, a feature which the inverse 
method cannot duplicate, figures 9 and 10 have been prepared. This 
calculation corresponds to the = 5 and y = 5/3 flow around the 
axisymmetric shape generated by rotating the principal meridian of the 
magnetopause about its axis and compares with an experimental test by 
Spreiter et al. (ref. 3). This particular profile shape contains a 
pronounced . dent with a concave corner in the vicinity of the neutral 
point. Spreiter et al. (ref. 3) have argued that the presence of a 
concave corner in the magnetopause is unrealistic and would not occur in 
.nature, but should be replaced by a free surface created by drawing a 
tangent line across the corner. In figure 9, we present results for both 
the concave and the free surfaces, denoted by solid and dashed lines, 
respectively. For the concave surface, a pocket of embedded subsonic flow 
is seen to form behind the embedded shock wave which was caused by the 
indented profile. If the surface is modified to the free surface shape, 
the embedded wave disappears completely as expected. The corresponding 
pressure coefficients on the magnetosphere boundary are presented in 
figure 10. For the concave surface, the shock wave is located on the body 
at approximately B = 80°, while for the free surface, the pressure 
coefficient displays a constant value as anticipated. Finally, we note 
that the calculation of a supersonic flow with an embedded shock and 
subsequent subsonic pocket provides a severe test for any blunt body 
procedure. The ability of the present code to provide convergent results 
for such a- flow demonstrates the capability for further extension and 
application to more generalized profiles than was heretofore possible. 

As a final indication of the range of solar wind flows over which the 
implicit code has been tested, we have prepared figure 11 which displays 
the variation of shock stand-off distance with oncoming Mach number and 
ratio of specific heats for three different magneto/ionopause shapes. 

These include the rotated equatorial trace of the magnetopause, and the 
ionopause shapes for H/R q = 0.01 and 0.5. Those shapes and the corre- 
sponding ranges of and y essentially span the entire range of interest 
of geometry and solar wind conditions for which the computer programs 
developed herein would be normally applied. 

In order to demonstrate the ability of the marching code to continue 
the blunt-body IMP solution downstream to some arbitrarily-specified 
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downstream location, and also to verify its stability and convergence for 
application to these classes of flows, we have calculated marching solu- 
tions for a variety of cases typical of solar wind interactions. Figure 

12 exhibits the location of the bow wave for flows with y = 5/3 and M = 

. 00 

3, 8, and 25 past the rotated equatorial trace of the magnetopause as the 
flow proceeds downstream from the nose region to x/D = -11. For this 
calculation, starting conditions for the marching code (SCT) have been 
provided by the blunt-body (IMP) code on the line x/D = 0.0, which is the 
usual location at which the two solutions are joined by the computer 
program. The marching code then determines the remainder of the flow 
field back to the specified downstream location. Since for the particular 
shapes considered herein the downstream flows are quite smooth, the 
marching calculation is very efficient (less than 30 sec., CDC 7600, OPT=2 
compiler) . Similar results are presented in figure 13 for M^ = 8 and y = 
5/3 flow past an ionopause shape with H/R q = 0.1. Those results have been 
carried downstream to x/R q = -20 since this distance is more typical of 
the region of interest for nonmagnetic planets. This is so because of the 
different scalings (D , R ) used for the magnetic and nonmagnetic applica- 
tions. For Venus and Mars, the location of the ionopause nose is at a 
distance only slightly greater than the planetary radius, as compared with 
the location of the magnetopause nose for the Earth which is of the order 
of 10 planetary radii. 

In order to illustrate the capability of the present procedures to 
determine streamlines, contour maps of various flow, properties , and 
magnetic field. lines- and contours, as well as to demonstrate the automated 
plotting capability for displaying these results, figures 14 and 15 have 
been prepared. Figure 14 illustrates the computer generated streamline 
locations and contour maps of density p/p^ and velocity ratio V/V^ for the 
complete near and far field flow about the equatorial trace of the magneto- 
sphere for M = 8 and y = 5/3. Based on this gasdynamic solution, figure 15 
exhibits the corresponding results for the magnetic field components 
(b/b^)^ and (B/B^)^ in the plane of symmetry. Both magnetic field line 
locations and contours are provided". In addition to demonstrating the 
overall smoothness of the computed results, these two figures illustrate 
the ability of the present techniques to provide the completely automated 
production of report quality plots of both gasdynamic and magnetic field 
properties for solar wind flows past axisymmetric magneto/ionopause shapes 
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CONCLUDING REMARKS 

i 

Theoretical analysis and associated development of computer programs 
were carried out for the purpose of developing computational techniques 
for predicting the interaction of solar wind flows with both magnetic and 
nonmagnetic terrestrial planets. Based on the identical theoretical model 
employed in previous work, i.e. the steady, dissipationless, magnetohydro- 
dynamic model for axisymmetric , supersonic, super -Alfvenic solar wind flow 
a new and more powerful solution procedure has been developed. The proce- 
dure is built upon an assemblage of computer codes, including: 

• Blunt-body code — to determine gasdynamic solution near obstacle 

nose 

• Marching code — to determine gasdynamic solution downstream of 

obstacle nose 

• Streamline code — to determine flow field streamlines 

• Contour code — to determine contour maps of flow and magnetic field 

properties 

• Magnetic field code — to determine frozen-in magnetic field 

• Plotting code — to plot selected flow and magnetic field results 

Comparisons are reported which demonstrate the accuracy of the present 
techniques by comparison with previously-established theoretical methods 
as well as limited experimental data. Furthermore, new results are pre- 
sented for a variety of solar wind flows which illustrate the flexibility 
and generality of the component methods. The computational procedures are 
fully automated and provide detailed flow field and magnetic field proper- 
ties in convenient outpuj; format, including an automatic plot capability. 
The programs are documented and are presented in a general user's manual 
included as part of this) report. 

With regard to improvements of the present techniques to provide a 
more accurate mathematical representation of solar wind flow, a number of 
items for further study are recommended. The first of these involves 
extension of the present capability from axisymmetric to three-dimensional 
nonaxi symmetric obstacle shapes representative of planetary magnetopauses. 
The second involves the merging of the present procedures for determining 
the magnetic field exterior to the magnetosphere boundary with a method, 
such as the Olsen-Pf itzer model, for calculating the confined magnetic 
field so as to provide a complete magnetic field representation. Lastly, 
the capability for determining numerical solutions of the complete magneto- 
hydrodynamic equations, rather than the frozen magnetic field approximation. 
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should be actively pursued. Although this extension represents a major 
advance in magnetospheric physics and computational fluid dynamics, in 
view of the success of the currently employed methods, its successful 
achievement appears highly promising. The component methods on which the 
present procedure is based are directly capable of extension in these as 
well as other directions. Moreover, the efficiency of the present calcu- 
lations indicates that the improved representation achieved by these 
improvements will not result in prohibitive computational costs. 


APPENDIX A 


A. 1 Introduction 

The purpose of this appendix is to describe the operation of the 
assemblage of computer codes which were developed in conjunction with the 
theoretical work presented in this report and organized into one program, 
and to provide sufficient detail to permit understanding and use of the 
program. The program computes the flow field of the solar wind about a 
terrestrial planet, using a procedure for the calculation of supersonic/ 
hypersonic flow about an axisymmetric blunt body. The corresponding 
frozen-in magnetic field is then calculated from the previously determined 
velocity and density fields. Streamlines and contour lines of various 
flow field properties and magnetic field components are also determined. 

A description of the general operating procedure of the program is 
given, with descriptions of input and output. The program is written in 
F0RTRAN IV and has been developed on a CDC 7600 computer. University 
Computing Company (UCC) Standard Plotting Software and Functional Software 
packages are used to produce automated plots. Files used, besides TAPE 5 
for INPUT and TAPE 6 for 0UTPUT, are TAPE1 for the plot file (system 
default) , TAPE 4 for input file for rerun option, and TAPE 9 for storing 
data for rerun. Typical run times for cases using the default parameters 
are 40 to 50 seconds, using the 0PT=2 compiler. For a case using the 
rerun option, which employs a previously calculated flow field, the run 
time is approximately 5 seconds. 


A. 2 Program Description 


For computational purposes, the flow is subdivided into two regions, 
as indicated in the sketch below, with the center of the planet as origin. 

iR 


Blunt Body 
Calculation 


Embedded 
subsonic region 

M >>1 


Marching Calculation 

-gBs *■ Bow Shock Wave 


use 



The region hear the nose of the magnetopause/ionopause includes all of the 
imbedded subsonic flow and part of the supersonic flow. An axisymmetric 
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implicit unsteady Euler equation solver is used to calculate this part of 
the flow field. Using the solution plane at x = 0 . 0 to provide starting 
conditions, the flow field in the purely supersonic downstream region is 
determined by integrating the steady Euler equations using a spatial 
marching procedure. Streamlines, the magnetic field, and contours are 
calculated using the entire flow field, distinguishing between the two 
regions as required by the different forms of the computational grids. A 
rerun capability is provided, where flow field data is read from a file 
written on a previous run, rather than repeating the blunt body and 
marching calculations. The computations proceed as shown in the sketch 
below, which provides an overall flow chart of the complete program. The 
program provides for several cases to be run consecutively. 
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A. 2.1 Calculation Procedure 

After reading in the number of cases in the run, each case is calcu- 
lated independently. Subroutine INPUT reads in all card input required 

for one case, viz. a title, flow conditions, obstacle geometry , calcula- 
\ ' 
tion and print control parameters, and desired contour values. The user 

* ‘V 

may supply the obstacle geometry in the form of a shape table for an 
axisymmetric body, or use one of the default shapes which are calculated 
internally by the program. These default shapes are the magnetopause 
equatorial trace and a constant scale-height .ionopause . for : which H/R q is 
specified by the user. The input is printed as the first item of output. 

A computational mesh in polar (R^,0) coordinates is established for 
the blunt body calculation; then, for the marching calculation, this is 
extended into a cylindrical (x,R) system, as indicated below: 



All lengths, x, R, R , are scaled so that the nose of the obstacle is at 

x = 1.0. For the default shapes, rays at equal angular increments of A0 

are used, starting at -A 0/2, up to 90° + A 0, where A0 = . 9 0 °/ (NBLUNT-1 . 5) 

and NBLUNT is an input parameter describing the number of angular mesh 

points to be used in the blunt-body calculation. Program default value is 

NBLUNT = 24, so that for the default mesh, A0 = 4°. The obstacle shape is 

determined by integrating the appropriate differential equation by a 

trapezoidal predictor-corrector method. For a user-supplied shape, the 0 

grid is determined by rays from the origin through the first NBLUNT points, 

and the reflection of the first ray in the x-axis. Values for R are 

P 

determined by dividing the line segments between the body and bow shock 
wave into NR-1 equal intervals. Thus, including the obstacle and bow 
shock wave, the grid forms NR arcs around the obstacle. A starting solu- 
tion for the blunt-body calculation is obtained by guessing a bow shock 
shape and by prescribing a Newtonian pressure distribution on the body. 
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Noting that the maximum entropy streamline wets the body, other flow 

properties on the body surface can then be calculated. An initial flow 

field is then established by linear interpolation between the obstacle 

and the guessed bow shock, where the Rankine-Hugoniot relations hold. The 

integration proceeds in time for ITER steps. The initial bow shock shape 

used for the magnetopause equatorial trace and for an ionopause with 

H/R q >_ 0.1 is a correlation shape depending on (M m , y, H/R ) and given by 

the parabola R = 6 , /6 -x/ZS - where 

p ioo 


6 

o 


1.0 + 1.1 


((y-l)M 2 

' CO 


+ ; 2)/(y+l)M 2 
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(0.9 + 0.5 H/R ) 
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i o 00 o 
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A 
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/ 
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x 0.78 


For a user-supplied obstacle shape and for an ionopause with H/R q < 0.1, 
the initial shock shape used is the curve R^ = Z(1 +A q (1 + 0.680* + 0 . 160 4 ) )'. 
Information on convergence, the final sonic line locations, and the body 
and final bow shock shape are printed from this calculation. 

The results at 0 = 90°, i.e. at x = 0.0, are used as starting condi- 
tions for the marching calculation, after proper normalization. For 
default geometries, the integration of the appropriate differential equa- 
tion is continued downstream at equal 0 increments to form a body shape 
table. The stepsize along the x-axis is recalculated at every IC0NST(49) 
with IC0NST(49) being set to 10. At each x-location, R^^ is determined 
by linear interpolation. The computational mesh is extended by adding the 
line perpendicular to the x-axis at each step, divided in the same manner 
as for the blunt nose. The calculation marches downstream with a maximum 
stepsize of 1.0 until the terminal location specified by the user has been 
passed. However, the number of steps is limited to 75, after which the 
calculation will end regardless of the x-location. The coordinates of the 
obstacle and bow shock are printed at each step. 


The grid coordinates and flow field values are written to a file, 
TAPE 9 , which may be saved to use as input for a later run.- This rerun 
option, which replaces constructing the computational mesh and performing 
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the blunt body and marching calculations with reading the rertfn input 
file, TAPE4 , is described in section A. 2. 2. 

The streamlines are calculated, in two, sections, following each of the 
flow field calculations. Using the, results of the blunt body calculation, 

i.e. the (x, R) grid coordinates, (R , 0) grid coordinates, density p/p , 

P "'ii °° 

and velocity components v x /v m and Vj^/v^ , the velocity magnitude | v | /v m 

and flow angle are calculated. Density p/p m and velocity Ivl/v^ are 

smoothed along the rays of constant-0, using a third degree least squares 

fit with respect to R . Streamlines are then calculated downstream to 

P 

x = 0.0, using the trajectory method, integrating through the velocity 
field by means of a third-order modified Euler integration procedure, 
using the grid locations on the bow shock as starting positions. The flow 
angle 4> = tan" 1 (v„/v v ) at each point is determined using bivariate linear 
interpolation in 0, then R '. Points for which 0 < 0° or 0 > 90° are 
discarded in the interpolation. 

The marching calculation provides (x, R) grid coordinates, and values 
of density P/P t , and velocity components v x / v t and v R /v t , where t denotes 
free-stream stagnation conditions. For compatibility with the blunt body 
solution, the flow field values are converted to p/p , v v /v , v D /v 
before calculating the resultant velocity magnitude Ivl/v^ and flow angle 
4) . The streamline calculation is continued downstream, employing the' same 
method as in the nose region. Starting positions on the shock wave for 
the streamline calculation in the marching- zone are set at equal R-incre-. 
ments, with a maximum of 50 streamlines calculated. The flow angle is 
determined using bivariate linear interpolation first in x, then in R. 

Along the symmetry axis, values of x, p/p^, and Ivl/v^ are determined 
by extrapolation, using a third-order Lagrangian polynomial in 0 on each 
arc of the computational grid. Exact values for the stagnation streamline 
are used where possible, viz. 

r • .... 

■ i ■ ’ i ■ 

at the bow shock 

P/Poo = (y+i)m^/ (y-DM^ + 2 j 

/ 

M/ v m = V(P/PJ 

\J 
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at the body surface 

rr f . ; v . ■ri 1/ (y~i) ] 

p/ Poo = (P/pJ sh o ck j [(Y+1 )M m J 2 /[4yM2 - 2 (y-1) J j 

Ivl/v^ =0.0 
x = 1.0 

Detailed flow field output may now be printed by subroutine FL0UT, with 
LPRFL as print control variable. In addition to grid coordinates, density, 
velocities and flow angle, values of temperature T/T^ and pressure P/P' 
are output, where 

T/T^ = 1 + (y- 1)/2]^1 - (Ivl/vJ 2 ] 

P / P 00 = (P/P„) ( T / T J : 

Streamline coordinates may also be printed by subroutine ST0UT , with LPRST 
as print control variable. 

The magnetic field is determined by separately calculating the 
components whose field lines are parallel, perpendicular, and normal to 
the flow, in the undisturbed solar wind/ These components are then added 
vectorially, the resultant being expressed in orthogonal (x,y,z) compo- 
nents. The angles in the free stream a p and a n between the magnetic field 
and the flow, as shown in the sketch below, are input variables. The 
components are calculated using the following formulae, in which e signi- 

'V 

fies a vector of magnitude e, in the direction of the component field line, 
and n the unit normal vector. 
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The magnetic field line vector component B„ which results from the inter- 
planetary component ,B w ^ that is parallel to the undisturbed solar flow has 
local magnitude given by (|v|/v ) (p/p- ), and the .same local direction <j> as 
the fluid flow. Determination of the normal magnetic field component 
requires calculation of R/R^, where is the free stream cylindrical : 
R-ordinate of the streamline through the point under consideration. This 
is calculated by linearly interpolating in the local radial cylindrical 
coordinate R between the streamlines ,' with R/R^ = 1.0 along the x-axis. 

The magnetic field vector component Bj. resulting from the interplanetary 
component which is perpendicular to the undistrubed solar wind flow 

requires the distance vector Ai/Ai^-, whose magnitude is A£/AJ l m and direc-< : 
tion is \p , where is the stretching factor of the perpendicular 

field at the point, and ip is the direction of the field line through the- 
point. The perpendicular field lines are determined by integrating /vdt 

T : . J . • * •• . ' ■ 'Xt 

along each streamline; using trapezoidal integration to locate points 
along the streamline at regular increments in time, At, starting at a 
perpendicular field line ahead of the bow shock. Values for Ai/Ai.^ and. ip . 
are calculated at the points where the perpendicular field lines, and 
streamlines intersect, interpolating only, along the field lines. A 
generalized quadrilateral.-interpolation scheme is then employed to 
determine hl/M^ and ip -at the computational grid points, using the quadri- 
lateral containing the point formed by the intersection of pairs of 
adjacent streamlines and perpendicular field lines. At the bow shock, an 
exact formula is -used, viz. ? 

(Ai/AS.^) 2 = 1 + cot 2 6(l+D 2 ), - 2D x csc0 x cot3 x cos(6-6) 
ip = 0 + sin -1 {D x cot0 x sin (0 — 6 ) / ( Ail/Ai^) } 
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where 

D 2 = 1 - 4 (M^sin 2 0-1) (YM^sin 2 0 + 1) / [ (y+1) 2 M^sin 2 0 j 
cot 6 = tan0 x | (y+l) m 2 / ^2 (M 2 sin 2 6 - 1)J - lj 
e . 

The values of at the grid points are smoothed using fifth-order 

least squares fit with respect to arc length along the arcs of the grid. 

The resultant magnetic field is expressed in orthogonal (x,y,z) components, 
where 


B /B = 
x 00 

cos 

a x 

n 

Jcosijl 

x COS 

me 

X 

ft 

a 

/B ) 

1 ' CO 

B /B = 

y' ~ 00 

cos 

a x 
n 

j^sint}) 

X cos 

CpXllBl 

/B ) 

' oo ' 

B /B = 

z 00 

sin 

a x 
n 

(B/BJ 

1 

n 




+ cosip x sin ci p x (IbI/B^) j.J 
+ sint|/ x sin a p x (| b j/B^) xj 


Magnetic field components may now be printed by subroutine B0UT, with LPRB 
as print control parameter. The magnetic field is not calculated when 
LPRB = .FALSE, and KBC0N=O . 

Contours are calculated for velocity \v\/v m , density p/p , and 
magnetic field components (B/B^),, and (B/B^Jx- The method used is a 
modified version of a procedure developed by R. Sorenson of NASA/Ames 
Research Center. The boundary is searched for intervals which bracket a 
contour point. Having found one point, the remainder of the contour is 
determined by 'walking' around the contour, searching at each step for the 
interval through which the contour line next passes, until a boundary 
point is reached. Then closed contours are found in a similar manner. 
Linear interpolation is used throughout the process. Note that since 
T/T m is a function of |v|/v ro only, velocity contours may also be consid- 
ered as temperature contours. The coordinates of the contour lines can be 
printed by subroutine C0N0UT, with LPRCON as print control parameter. The 
section of the program which creates the plots is accessed only when the 
control parameter LPLOT = .TRUE. 


42 



APPENDIX A 


A. 2. 2 Rerun Option 

The rerun option is used when LRERUN = .TRUE. The blunt body and 
marching calculations are replaced with reading the grid coordinates and 
flow field values from the rerun file, TAPE 4 , which contains data written 
to TAPE9, then saved, on a previous run. Different values for any param- 
eter not used in the flow field calculations may be specified, viz. 
contour values, plot length, magnetic field angles, and output options. 
Values of AMACH, GAMMA, and HR0 are required input, to ensure that the 
input rerun file does contain the case desired for rerun. If the geometry 
is user-supplied, the body shape table will be read from TAPE 4 , and should 
not be input from cards. 

After reading the card input, AMACH, GAMMA, and HR0 are tested 
against values from TAPE 4 . The grid coordinates and flow field values 
from the blunt body calculation are read in, then smoothed, and stream- 
lines calculated for this region, as previously described. The results 
of the marching calculation are then read, and the streamline calculation 
continued downstream. The calculations then proceed as described in 
section A. 2.1. 

A run must not contain more than one case which uses the rerun option 

A. 2. 3 Program Limitations and Precautions 

The program makes some assumptions about the geometry of the obstacle 
shape around which flow is to be calculated, and about the flow field. 

The obstacle shape is assumed to be monotonically increasing in R, going 
downstream. The nose of the obstacle is at x = 1.0. The origin of the 
(x, R) coordinate system is the center of the planet. Obstacle shapes 
with sharp corners should be avoided. The flow at x = 0.0, which is used 
as the starting plane for the marching calculation, must be such that the 
axial (x-direction) Mach number is purely supersonic. This is a basic 
requirement of the marching code and is associated with the manner in 
which the code advances the solution. That condition is violated for 
ionopause shapes with H/R0 > 0.5, so that with the present joining proce- 
dure employed (i.e. at x = 0.0) those shapes cannot be treated. In the 
magnetic field calculation, the first streamline is assumed to be inside 
the arc described by the grid points immediately off the body, downstream 
of x = 0.0. To reduce computational costs, a grid using NR = 10 may be 
used, in which case a lower .value of CN may be required. This would 
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reduce the running line by approximately 40 percent. A free stream Mach 
number less than 3.0 is not advised. 

A. 2. 4 Convergence Criteria for Blunt Body Calculation 

The output provides two measures of the convergence of the blunt body 

calculation. The RMS of shock speed and maximum shock speed are printed 

at each iteration. These quantities should both tend to zero. A value 

for RMS of shock speed, of ; 

^RMS 

q RMS < ^ * 10 " 3 

where y ‘is the specific heat ratio, arid is the free stream Mach number, 
usually indicates a converged solution. The RMS of error in enthalpy, HT , 
should be less than 1 percent, with the maximum enthalpy error also of that 
order . 

The Courant number, CN, determines the time step size used by the 

v * ’ f 

calculation. A value not greater than the default of 3.0 should be used. 
For low Mach numbers or a coarser mesh than the default, a lower value may 
be preferable. If the default value does not generate a converged solu- 
tion, or if the error message from subroutine SH0CK is printed, try 
lowering CN in increments of 0.5. to find a better value of CN. User- 
supplied bodies may also require a lower Courant number. 

A .3 Description of Input 

This section describes the card input for the program. An alpha- 
betized dictionary of input variables is provided, defining the variables, 
listing default values and limitations. A discussion of the preparation 
of the card input is then presented, followed by a description of the 
input card format. 

A. 3.1 Dictionary of Input Variables 

AMACH free-stream Mach number; 3 . 0 <_ AMACH £ 25.0 is recommended 

ANGN the angle, in degrees, measuring the deviation of the free- 

stream magnetic field from the plane in which B „ and B , lie; 
equal to tan^B^//^. * | +1^1*) 

ANGP the angle, in degrees, measuring the deviation of the in-plane 

magnetic component (B + B ) from the direction of flow; equal 
to tan - 1 (B /B ) ^°°" ' XjC ° 1 

' coj/ oo i r 
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BC0N (I) 
CN 

GAMMA 

HR0 

ITER 

KBC0N 

KRC0N 

KVC0N 

LPLOT 

LPRB 
LPRC0N 
LPRFL 
LPRST - 


KBC0N-dimensional array specifying values to be used for mag- 
netic field strength contours 

Courant number used for blunt body calculation; program default 
value is 3.0 

free-stream ratio of specific heats 

obstacle geometry indicator: 

HR0 > 'O'— ionopause '• with H/R Q = HR0 
HR0 = 0 — magnetopause equatorial trace 
HR0 < 0 — geometry is user-supplied 

integer, number of iterations for blunt body calculation; 
program default value is 300 

integer ,. number of values specified for magnetic field contours; 

0 £ KBC0N <_ 20 

■. ’ 1 7 . • ‘ I . K • 

, integer, number of values ^specified for density contours; 

0 < KRC0N <_ 20 

integer, number of values specified for velocity contours; 

0 £ KVC0N <.20 

logical variable indicating whether to create plots or plot file 

FALSE .—no 
TRUE -‘yes 

logical variable indicating whether to print magnetic field „ 
output 

FALSE — no 
TRUE — yes 

logical variable indicating whether to print coordinates of 
: contour lines - ' 

FALSE - no 
■ TRUE — yes 

logical variable indicating whether to print detailed flow field 
output 

FALSE — no 
TRUE - yes 

logical- variable indicating whether to print coordinates of 
streamlines 

FALSE — no- __ .i 

TRUE — yes 


45 



APPENDIX A 


L RE RUN 

NBLUNT 

NB0D 

NCASE 

NR 

"RC0N ( I ) 
RR (I) 
TITLE' 

VC0N ( I ) 

XCALC 

XPL0T 

XX(I) 


logical variable indicating whether this case uses rerun option 

FALSE — perform blunt body ;and marching calculations 
TRUE — read results of a previous calculation from TAPE 4 

integer, number of angular mesh points for blunt body calcula- 
tion; for user-supplied geometry, XX (NBLUNT-1) =0 . 0 ; program 
default value, and maximum; is 24 

integer, number of points in body shape ' table when geometry is 
user-supplied; 1 <_ NB0D £ 100 

integer, number of cases to be run consecutively; NCASE >_ 1 

integer, number of radial mesh points; program default value, 
and maximum, is 19 

KRC0N — dimensional array specifying values to be used for 
density contours ‘ 

j ■ . ■ 

NB0D — dimensional array representing the R-locations, in 
cylindrical (x,R) coordinates, of the user-supplied body shape 

descriptive heading of the case, to be printed on the first 
page of output; may contain ’up to 80 characters, including 
blanks 

KVC0N — dimensional array specifying values to be used for 
velocity contours 

terminal downstream x-location for marching calculation of flow 
field; XCALC < 0.0; program default value is -1.0 

terminal downstream x-location for calculation of streamlines, 
magnetic field, and contours; XCALC <_ XPL0T < 0.0; program 
default value is -1.0 

NB0D — dimensional array representing the x-locations, in 
cylindrical (x,R) coordinates, of the user-supplied body shape 


A. 3. 2 Preparation of Input Data 

The card input for a run consists of one card containing the number 
of cases to be run consecutively. Item 0, followed by a set of input for 
each case, Item 1 through Item 7, and Item 8 if required. Where a default 
value is to be used, the input field should be left blank. For each case, 
all required variables which do not assume their default values should be 
specified. The input format for. all cards is described in section A. 3. 3. 

Item 0 — This item consists of one card, containing the number of 
cases in this run, NCASE. 

Item 1 — This card provides' identification of the case, TITLE, which 
is printed on the first page of the output for this case. 
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Item 2 — This card contains information on the flow conditions and 
body geometry, and parameters required for the blunt body and marching 
calculations. AMACH, GAMMA, and HR0 must be specified for each case. For 

the. rerun option, the values are tested against the values from the rerun 

file. The parameters XCALC, NR, NBLUNT, CN , ITER are used only when the 

flow field is to be calculated. These variables each assume a default 
value if the input field is blank. 

Item 3 — This item consists of one card containing the rerun indi- 
cator, LRERUN, and the output control variables LPRFL , LPRST , LPRC0N , 

LPRB , and LPL0T . 

Item 4 — This card contains the variables XPL0T , ANGP, and ANGN. The 

value for XPL0T is changed by the program to be the x-location of the 

marching calculation immediately upstream of the input value for XPL0T . 

The angles describing the deviation of the magnetic field from the flow, 

ANGP and ANGN, are not required when LPRB=. FALSE, and KBC0N=O , since the 

magnetic field is not calculated under these conditions. ANGP is the 

angle between the vectors (B . + B „ ) and v , while ANGN is the angle 

between B and (B + B ) , where B , B , B^ are the components of the 
0 ,°° \,°°->- 'v 00 " v»i« 'v»a v»n r 

free-stream magnetic field, B , which are parallel, perpendicular, and 

Oi 00 

normal to v , and are as indicated in figure 4 of the text. The two 
angles ANGP and ANGN fully determine the half plane for which the magnetic 
field is to be calculated. The magnetic field for the other half of the 
plane may be calculated by rerunning with the sign of ANGP reversed. When 
(B + B ) = 0, ANGN = ±90°, ANGP = 0°; and, when B =0, ANGN = 0°. 

\°°JL ^0011 ry/jon 

Item 5 — This item contains the values for the velocity contours. 

The first card contains KVC0N, the number of values specified for VC0N. 

If KVC0N=O, only this card is required, and no velocity contours are 
calculated. If KVC0N > 0, the contour values are then read. Up to three 
cards may be required to accommodate the values, eight per card, maximum 
of 20. The contour values should be monotonically increasing, with at 
least one value within the range of the magnitude of the velocity in the 
region for which contours are to be calculated. 

Item 6 — This item contains the values for the density contours. The 
description is similar to that for item 5, with KRC0N being the number of 
values specified, and RC0N the array of values. 

Item 7 — This item contains the values for the magnetic field 
contours. The description is similar to that for item 5, with KBC0N 
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being the number of values specified, and BC0N the array of values. Note 
that the same contour values are used for the parallel and perpendicular 
components . 

Item 8 — This optional item is required when HR0 < 0.0 and LRERUN= 
.FALSE.’, and contains the body shape table for the user-supplied geometry. 
The first card contains NB0D, the number of points in the shape table. 

The next NB0D cards contain the cylindrical (x,R) coordinates of these 
points,^ (XX ( I ) , RR ( I ) ) , one point per card. The points supplied by the 
user determine the 0-spacing of the mesh used for the! blunt body calcula- 
tion. The first point should be near,, but not on, the x-axis. A suggested 
location is such that the 0-spacing between the first point and the x-axis 
is half the 0-spacing between the first, two points. The blunt body calcu- 
lation adds a point which is the reflection about the x-axis of the first 
point in the body shape table. The (NBLJJNT-1 ) point should be at x=0.0. 
The BLUNT 1 " point is also used to create the grid for the blunt body 
calculation. ; The coordinates must be normalized so that the planet center 
is at ( : 0.,0.) and the nose of the body at; (l.,0.). 

-A. 3. 3 Format of Input Data . \ 

Four format types are used for the input data. For real numbers 
(F-format), a decimal point is required ., Integers (I-format) should be 
right-adjusted ^in the field,. For logical variables (L-format) , the first 
non-blank character in the field,; which should be ’T’ or ’F’, determines 
the value. Note that a blank input field is interpreted as ’FALSE'. The 
title, which, is; in A-f brmat,,. ; may jcontain any valid character. 

A description of the card format of the input data follows, with 
item numbers corresponding -to those i : n . section A. 3. 2: 
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Input Data Card Fbrmat (Continued) 
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A. 4 Description of Output 

This section describes the output of the computer program. The 
contents of each output item are specified and discussed. The printed 
output consists of six items, four of which are optional, with input 
control parameters. Plotted output is also optional. 

The first output item consists of the input data. The title for the 
case is printed first, followed by all input variables, with identification 
of each variable. Default values are printed as if they were input. 
Parameters CN, NR, NBLUNT , ITER for the blunt body calculation and XCALC , 
the terminal location for the marching calculation, are printed only when 
the flow field is to be calculated. When the geometry is user-supplied, 
the input body shape table is printed. For a default geometry, the body 
shape is indicated by the description "default magnetopause coordinates — 
equatorial trace" or "default ionopause coordinates — H/R0 = ". 

The second output item is not printed when LRERUN= . TRUE . From the 
blunt body calculation, the shock speed at each iteration, the final 
enthalpy error, final sonic line location, and body and final bow shock 
shape are printed. For convergence criteria of this calculation, the 
downstream location and body and shock ordinates are output. There is no 
control variable allowing the user to suppress this item of output when 
the flow field is calculated. 


Detailed flow field output is the third item, and is printed only 
when LPRFL= . TRUE . Coordinates are labeled as X/D, R/D, RP/D, or X/RO, 
R/RO, RP/RO, to emphasize that distances are normalized by the distance 
from the center of the planet to the nose of the body, D for the magneto- 
pause, RO for an ionopause. Along the symmetry axis, the values printed 
are velocity magnitude V/VINF, density RH0/RH0INF, temperature T/TINF, and 
pressure P/P INF. Over the rest of the flow field values are also given 
for velocity components VX/VINF, VR/VINF , and flow angle <j> . Note that the 
flow angle is the deviation of the flow about the obstacle, and so 
0° < <j> < 90°. 


The next output item is the (x,R) coordinates of the streamlines. 

For the blunt body region, the (Rp,0) coordinates of the starting position 
on the bow shock wave are also given. This item is printed only when 
LPRST= . TRUE . 




/ 
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The magnetic field components are then printed,, if LPRB= .TRUE . . The : . 
location of each point is defined in (R p ,0) coordinates for the blunt body 
region, and (x,R) coordinates for the downstream marching region. The 
components along field lines- parallel, perpendicular, and normal to the 
flow in the free stream are printed as B/BINF (PARALLEL) , B/BINF (PERP) , 
B/BINF (N0RMAL) . The orthogonal (x,y,z) components of the resultant are 
printed as BX/B INF (RESULTANT) , BY/BINF (RESULTANT) , BZ/BINF (RESULTANT) . 

The magnetic field in the symmetry (x,y) plane, defined by the vector sum 
{(B/B oo ) i( .+ (B/B^) }, is also, printed, and is given by the magnitude 
B/BINF (IN-PLANE) and direction. B-ANGLE ( IN-PLANE) of the vector. The 
B-ANGLE is- the angle between the magnetic field and the flow, and 
generally has the same sign as ANGP . 

The last item printed is the (x,R). coordinates of the contours, for 
which LPRC0N is the logical control variable. Noting that temperature and 
velocity contours coincide, the corresponding value of T/TINF is printed 
along with V/VINF for the velocity contours. There are three nonfatal 
error messages which may occur — see section A. 5. 

The program also has the capability of producing plotted output, 
using UCC plot routines AXIS, CHAR,, DASH, D0TLN, ENPLT, GREEK, MATH, 
NUMPLT, PL0T, PLTLN , P0LAR, SCALF , VECT0R. Plots are not produced when 
LPL0T= . FALSE . The plots provide pictorial representation of the stream- 
lines and contours, with a maximum of six frames produced. The first 
frame is a plot of the streamlines. Then- contour plots of velocity, 
temperature, and density are drawn. The last two frames are contour plots 
of the parallel and perpendicular magnetic' field components , with field 
lines added to indicate direction. Frames' are omitted when no contour 
levels are provided.. 

A. 5 Program Error Messages 

This section lists the messages printed by the program, 
what action should be taken by the user. 


and indicates 



APPENDIX A 

( 1 ) 


FR0M CARDS 
FR0M TAPE 4 

The first three parameters of item 2 of the input for a case using 
the rerun option should agree with those used when creating the file. The 
tolerance used in comparing the values is 10 -5 . For a user-supplied 
geometry, it is sufficient for both 1 values of H/R to be negative.,' 

. ... ; ■ ! I . ° . . = • ' - ' . ' ' ; 

(2) ***** EXECUTI0N 1 TERMINATED ***** 

ARRAY 0F C0NT0UR VALUES IMPR0PERLY SPECIFIED ‘ 

When specified, the contour values should be -monotonicallyincf easing, 
with at least one value in the range of the velocity; density, or magnetic 
field strength for the region under consideration: This' error does not 

inhibit generation of the rerun file. • ' •• * 

(3) C0NT0UR SEARCH AB0RTED. TABLE 0VERFL0W IN NAD , 

The program allows for 29 contour’ lines’ to r be found, storing the 

# . .. V • '*.^P 

starting address of each contour in: array NAD. This message’ indicates 
that’ at least one more contour line could be found. If the user requires 
all the contours of the levels specif ied , the case should be rerun in two 
parts. Otherwise, reduce the number of contour levels specified. 

(4) C0NT0UR SEARCH AB0RTED — TABLE- 0VERFL0W IN. (X,Y) 

The contour lines may be described by up to 1000 points, stored in 
arrays X and Y. This message indicates that more points would be required 
for the contour lines requested. The last contour line found will be 
incomplete. As with (3), either reduce the number of contour levels or 
run as two cases. ; 

(5) NEGATIVE PRESSURE DETECTED BY SH0CK AT J= 

PN= P0= PTAU= 

This message is printed when a negative pressure has been calculated 
at the shock on this iteration, at radial location J. The quantities 
printed are: PN, the pressure calculated on this step; P0, the pressure 

from the previous step; and PTAU, the partial derivative of pressure with 


• ***** EXECUTI0N TERMINATED 1 ***** . 
RERUN DATA 0N TAPE 4 . D0ES N0T AGREE. 
WITH CASE SPECIFIED 0N ’CARD INPUT: 
MACH N0 . GAMMA - : M/RO 
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respect to time. This condition indicates that the shock wave motion is 
too extreme. Lowering the value of CN, and thus reducing the time step., 
may remove the problem. 

The following messages (6) -(10) usually result from using an obstacle 
geometry which is in some way too severe for the program to handle in its 
present form. The obstacle slope may be sufficiently high at x = 0.0 that 
the axial Mach number becomes subsonic in the starting solution for the 
marching calculation, or there may be a sharp corner in the profile. 

Check input, particularly free stream Mach number and body geometry. 

(6) NEGATIVE PRESSURE 0N B0DY DETECTED BY BNDRY , PB= AT J= 

This message indicates that a negative pressure on the body, PB, has 
been calculated at radial location J. 

(7) NEGATIVE PRESSURE 0R DENSITY 0N B0DY DETECTED BY BNDRYM AT X= 

PB= RH0B= VXB= VRB= 

The program makes internal corrections when this condition occurs, 
resulting pressure PB , density RH0B, and velocity components VXB and VRB . 

(8) NEGATIVE SIGMA-BAR-1 IN EIGENM INDICATES SUBS0NIC FL0W AT 1= 

(9) NEGATIVE SIGMA-BAR-2 IN EIGENM INDICATES SUBS0NIC FL0W AT 1= 

These messages are printed when subsonic flow is detected by the 
marching calculation. The computed stepsize for this region will be quite 
small . 

i 

(10) BODY TURN STOPPED AT M2=100 

This message indicates that the body has a sharp corner, which has 
been limited to. 100° when being transformed. 

A. 6 Sample Case 

A sample case is presented in this section. Figures A.l, A. 2, and 
A. 3 illustrate the input data, portions of the printed output, and the 
plotted output, respectively. 

The sample case is run alone and is set up so as to produce all 
possible output. Flow is to be calculated about the rotated equatorial 
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trace of the magnetopause to a downstream location of x = -10.0, with 

= 8, y = 5/3. • Streamlines, magnetic field components, and contours are 

desired to a downstream location of x = -5.0. The magnetic field forms 

angles of a = -45° , a = 45° with the solar wind flow. Contour values 
r p n 

are specified for all quantities. 

The input data is tabulated in figure A.l, with item numbers corre- 
sponding to those in sections A. 3. 2'- and A. 3. 3. The first card, item 0, 
indicates that there is one case to be run. The remaining 13 cards provide 
the data required for this case. Item 1 contains the identifying title. 

On the next card, item 2, values are specified for AMACH, GAM, HR0, and 
XCALC . The other data fields are left blank to indicate that the default 
values will be used. The values of the logical variables of item 3 specify 
that the flow field is to be calculated, and that full printed and plotted 
output is^to be produced. Item 4 defines the plot length to be -5.0, and 
the magnetic field angles as ANGP = -45.0 and ANGN = 45.0. Items 5 , 6, 7 
specify the contour levels to be used — 12 for velocity and temperature, 

10 for density,, and 12 for magnetic field strength. Item 8 is omitted, 
since the obstacle geometry is one of the default shapes for which the 
coordinates are calculated internally by the program. 

Figure A. 2 presents portions of the printed output from the sample 
case using the data deck of figure A.l. The full printed output for this 
case is approximately 5,000 lines. The first item of printed output, 
figure A. 2(a), is the input variables. Then in figure A. 2(b), the output 
from the blunt body calculation is shown, omitting shock speed values for 
iterations 6 to 295. Note that the convergence criteria described in 
section A. 2. 4 are satisfied. Figure A. 2(c) provides the output from the 
marching calculation. Figures A. 2 (d) , (e) , (f ) show the detailed flow fi^eld 
output in four areas: near the symmetry axis; near x = 0.0, when the two 

calculation schemes are joined; around x = -5.0, the value of XPL0T ; and 
the last two locations calculated. An example of the streamline output is 
shown in figures A.2(g),(h). Figures A. 2 (h) , (i) , ( j ) provide the magnetic 
field components at the same locations chosen to illustrate the flow field 
output. Examples of the printed output of the contour lines are shown in 
figures A. 2 (k) , ( 1) , (m) , choosing two contour lines for each quantity. Note 
that for ( I Bl/B^) , the different contour lines for the same contour level 
are considered to be two separate and distinct lines. 

Figure A. 3 shows the six plots which are produced by the program for 
this case. 
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Figure A-l. - Card input for sample case. 
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• 

• 7U 


• 8CC 
3. SCC 


*500 

2.K* 


FIGURE A. 2 


abbreviated output for sample case 


sample case default magnetopause equatorial trace 

INPUT VARIABLES 

FRftSThEAH MACH NG. • 8,0u 

SPECIFIC He AT RATIO • 1.667 

OBSTACLE GEOMETRY! DEFAULT MAGNETOPAUSE COORDINATES - EQUATORIAL TRACE 

PARAMETERS FOR HUNT 30DY CALCULATION 

NO. OF RADIAL HE S N POINTS - i? 

NO. OF ANGULAR HE SH POINTS • 2 A 

COL'kAnT NUMBER - 3.#£- 

NO. OF ITERATIONS • Hi, 

TERMINAL DOWNSTREAM LOCATION FOR MARCHING CALCULATION — H.CC 

TERMINAL DOWNSTREAM LOCATION FQ« PLOTTING • -S.wJ 

IN-PLANc DEVIATION OF MAGNETIC FIELD FROM FLOW • -AS. 00 DEGREES 

DEVIATION JF MAGNETIC FIELD FROM Plax€ If FLO* • AS. 00 DEGREES 

LRCRLN • F 
IPftFL « T 
LPRST - T 
IPRCON • T 
LPRB ■ T 
LPLOT • T 


VALUES SPECIFIED FOR CONTOUR CALCULATION 


12 CONTOUR LEVELS FOR VELOCITTi 

• 2s » .30.. . AOO 

*7jJ «tW *o5** 


:c CONTOUR L c WEL S FOR DtftSITM 

1.IC9 4.65: 2.L33 

3.MC A.l!C- A.2Q1 


• SOD .639 

• daj .9 >5 


2. SOS 3.000 


It CONTOUR LEVELS FOR MAGNETIC FlUO STREN6T H » 

• 75C I. COO 1.250 2.503 2.000 

3.CCL 3.SC3 A* 1 03 5. Cat 6.030 


Figure A. 2(a) 
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dLUit 30DT CUCJILTIO* 


iteration 

1 

RRS 

Of 

SHOCK 

SPEED* 

1.2458 c^ : *Z 

RAXIHUR 

SHOCK 

SPEEO- 3.6519E-C2 

AT 

J-26 

ITERATION 

2 

RRS 

Of 

SHOCK 

SPEED- 

1.2445E-D2 

RAXlHUR 

SHOCK 

SPEED- 3.1746E-02 

AT 

J-25 

ITERATION 

3 

RRS 

OF 

SHOCK 

SPEED- 

2.1994E-02 

raxIrjr 

SHOCK 

SPEED- 3.2075E-02 

AT 

J-25 

ITERATION 

4 

RRS 

Of 

SHOCK 

SPEED- 

3.4842E-D2 

Raxihur 

SHOCK 

SPEEO— 5.037CE-C2 

AT 

J-l 6 

ITERATION 

5 

RRS 

r*f 

SHOCK 

SPEED- 

4.94442-52 

maximum 

SHOCK 

SPEED—7.2259E-02 

AT 

J-26 

ITERATION 

296 

RRS 

Of 

SHOCK 

SPEED- 

3. 79 JO t-Ji 

RAX I HUH 

SHOCK 

SPfccU* 9.6232fc-13 

A 1 

J-l? 

ITERATION 

297 

RMS 

Of 

SHOCK 

SPEED- 

3.7524E-D3 

RAX1RJH 

SHOCK 

SPEED- 9.U11E-C3 

AT 

J-15 

ITERATION 

299 

RrS 

Of 

SHOCK 

SPEED- 

3.714GE-03 

HAXIRUR 

SHOCK 

SPEEO— 9.1274E-03 

AT 

J*17 

ITERATION 

299 

RRS 

Of 

SHOCK 

SPEED- 

3a6751c*>3 

RAX 1 HUH 

SHOCK 

SPEED— 8.9772E-D3 

AT 

J-17 

iteration 

30C 

RMS 

Of 

SHOCK 

SPEcD- 

3.6352E-D3 

HAXIHUH 

SHOCK 

SPEED— 8.6374E-D3 

AT 

J-14 


FINAL 

SONIC 

LINE LOCATION 

XSL- 

*5952 

RSL- 

.9999 

XSL- 

.617: 

RSL- 

.95 71 

XSL - 

.6297 

RSL - 

• 97 C 2 

XSL- 

• 654« 

RSL- 

.9812 

XSl- 

• 6792 

9 SL - 

• 9929 

XSL- 

• 6959 

RSL - 

.9985 

XSL- 

.7187 

RSL* 

1.CD93 

XSL- 

.7423 

RSL- 

.•.•i': e; 

XSL- 

• 7(65 

RSL- 

:.*i n 

XSL ■ 

.7919 

RSL- 

1.MC4 

XSL- 

.8172 

RSL- 

1.0D91 

XSL- 

. 6«»3? 

RSL- 

J.x\ c2 

XSL- 

• 3fc 99 

RSL- 

l.Cw24 

XSL- 

.9967 

RSL- 

.99 51 

XSL- 

.9290 

RSL- 

.9871 

XSL- 

• 951b 

RSL- 

• 977l 

XSL- 

.9795 

RSL- 

. 96t 1 

X5L- 

1.0177 

RSl- 

.95 29 

XSL- 

1.037C 

RSL ■ 

.93 7? 


«« X ee*Of IM HT * ,S822t»C3 *HS Of X EURO* IN NT ■ .M7H-31 >»» 


BOOT AND FINAL SHOCK S.NAft 


J 

X(BOOY) 

R (BOOT) 

X(SHOCK) 

R (SHOCK ) 

2 

-.9995 

• C 399 

-1.27)3 

• J4*9 

3 

-.9957 

• 1 C9 6 

-1.2662 

• 1331 

9 

-.9879 

.1792 

-1.2532 

• 2*ib 

5 

-.9769 

.2939 

-1.2956 

• 3136 

6 

-.9 tlC. 

• 1122 

—1.2 232 

.3994 

7 

-.9917 

.28» 5 

-1.2C82 

• 4882 

8 

-.918* 

• 99dL 

-1.1839 

• 5772 

9 

-.8915 

.519? 

-1.1537 

. 666. 

Xu 

-.8607 

• ;et 5 

—2 *22 st 

.7554 

11 

-.826C 

.(953 

—1.0619 

.8449 

12 

-.7879 

• 7l b9 

-j .C383 

.9349 

13 

-.7996 

.7713 

-.9898 

1.0250 

19 

-.6939 

.(323 

-.9167 

1.1163 

15 

-•( 48C 

• b 9.9 

-.8777 

1.2282 

16 

-.5935 

• 9999 

-.8130 

1.K-1J 

17 

-.535r 

i.rc (2 

-.791b 

1.3951 

18 

-.9723 

i.((C9 

-•6t*r 

1.4924 

19 

— .9059 

2.1137 

-.5783 

1.5886 

2D 

— • 3?9C 

2 »It*i8 

-.9«9'9 

1.6892 

21 

-.253. 

1.2238 

-.33*9 

1.7922 

22 

-.1772 

2 *2 611 

—.2668 

1.3937 

23 

— -091 3 

i.SLtf 

-.19;5 

2.JT88 

29 

. A 00C 

1.3991 

-•©Dit- 

2.1245 


Figure A. 2(b) - continued 
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fU*CHlNG CALCULATION 


* NO, 

OOUNSTRE Aft LOCATION 

BOOT ORDINATE 

SHOCK ORDINATE 

1 

-•L3C2 

1.2622 

2.1483 

2 

-.0604 

1.3751 

2.1719 

3 

9( 6 

l. 3374 

2.1951 

A 

l?Cb 

1.3993 

2.2181 

5 

-.1510 

1.41C8 

2.2406 

6 

-•1812 

1.4219 

2.2633 

7 

-.2114 

Z.4326 

2.2855 

8 

—,241b 

1.4430 

2.3074 

9 

-.2716 

1.453J 

2.3292 

10 

-.3253 

1.4699 

2.3671 

11 

-.3786 

1.4359 

2.4043 

12 

-.4322 

4. 5'.' 8 

Z.44V9 

13 

-.4ei? 

1.5149 

2.4769 

14 

-.5292 

1.5292 

2.5123 

15 

-.5926 

1.5418 

2.5472 

16 

-.64fcj 

1.5526 

2.5815 

17 

-.6996 

1.5&38 

2.6153 

18 

-.7530 

1.5743 

2.6487 

19 

-.8065 

1.5843 

2.6816 

2D 

-.6964 

1.5990 

2. 7359 

21 

-.9863 

1.6139 

2.7891 

22 

-l.f 7fci 

1.6267 

2.8412 

23 

-1.1661 

1.63t4 

2.8922 

24 

-1.2560 

1, 6491 

2.4423 

25 

-1. 3459 

1.6587 

2.9914 

26 

-1.4358 

1. 6676 

3**397 

27 

-1.5257 

1.6757 

3.0672 

28 

-1. tlSt 

1.6331 

3.1339 

29 

-1.7' 55 

1.6899 

3.1799 

30 

-1.6527 

1.6999 

3.2537 

31 

-1.9999 

im 7086 

3.3257 

32 

-2. 2471 

1.7161 

3.3962 

33 

-2.2943 

1.7228 

3.4651 

34 

-2.4415 

1.7287 

3. 5327 

35 

-2.5887 

1.7338 

3.5989 

36 

-2.7359 

1.7384 

3.6640 

37 

-2.8621 

1.7425 

3.7278 

38 

-3.C3C3 

1.7462 

3.790b 

39 

-3.1775 

1.7494 

3.8524 

4? 

-3.4261 

i.7542 

3.9545 

41 

-3.6747 

1.7582 

4.0541 

42 

-3.9232 

-4.2718 

1.761b 

4.1513 

43 

1.7545 

4.2465 

44 

-4.4ZC4 

1.7670 

4,3396 

43 

-4.6689 

l.ToOl 

4.93 l«r 

46 

-4.9175 

1.7711 

4.52J6 

47 

-5.1661 

1.7726 

4.6086 

48 

-5.4147 

1.7741 

4.6951 

49 

-3.t632 

1.7753 

4.7802 

50 

-6.C798 

1.7771 

4.9198 

51" 

-6.4963 

1. 77b6 

5.0561 

52 

-6.912b 

1.7799 

5.1894 

53 

-7.3293 

1.7809 

5.3190 

34 

-7.7458 

1.7918 

5.4476 

55 

-8.-1623 * 

1.7026 

5.5730 

36 

-6.*70S 

1.7833 

5.6961 

57 

-8.9434 

1.7839 

5.8173 

58 

-9.4U9 

1.7944 

5.9364 

59 

-9.8284 

4. 7648 

6. 6*538 

60 

-10.5CCC 

1.7054 

6.2394 


Pigure A. 2(c) - continued 
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DETAILED FlQ* FIELD OUTPUT 


FLOM 

FIELD VALUES 

EXTRAPOLATED 

TO SYNHSTRY 

AXIS# THETA 

- :.00 DEGREES 

I 

X/0 

V/VINF 

RHO/RHOINF 

T/TINF 

P/PINF 

1 

2 * (• Va> C 

O.Uf i- 

4,2291 

22*3333 

94.4497 

2 

1.0150 

• Gle4 

4.2200 

22.3276 

94.2214 

3 ' 

l.f.3 Ji 

• t 3i 7 

4.2156 

22.3133 

94.0645 

4 

i. 

• C 451 

4.27 o9 

22,2899 

93.8164 

5 

1.0602 

•C547 

4.2998 

22.2572 

93.476i 

& 

1.C752 

• C745 

4*1883 

22.2149 

93.043C 

7 

1 .t 903 

• C894 

4.2744 

22.1629 

92.5177 

8 

1.1C53 . 

.1043 

4.1562 

22.1022 

91.9009 

<9 

I.! 2?4 

.2192 

4.1395 

22.0300 

91.1941 

10 

1 .2 35 A 

.2341 

4.2 1 85 

22. 9494 

90.3992 

11 

1.1515 

. 1 49C 

4.0951 

21.3599 

39.5189 

12 

1 .1656 

..637 

4.C 093 

22.7518 

33.5562 

13 

2.1 60 5 

• 276 2 

V.C422 

22.5558 

87.5143 

24 

1 . 1 95 fc 

• 1925 

4.0106 

21.5424 

36.3989 

15 

2.2126 

• 2(66 

3.9777 

21.4225 

95.2129 

16 

1.2257 

• 22(4 

3.9424 

21.2969 

83.96*8 

17 

1.2AP7 

.2339 

3.9048 

21.1666 

82.6508 

! 8 

l.»55fe 

.2469 

3.6647 

22.0326 

81.2854 

19 

1.2706 

.2627 

3.8209 

20. 8721 

79.7500 


. r 




FLO* FIELD 

VALUES FRON 

3LUNT BODY CALCULATION 




ANGULAR 

LOCATION NO. 

2# AT THETA 

• 2.> GC 

> OEGREES 







I 

RP/D 

R/D 

X/D 

VR/VINF 

VX/VISF 

FLOW ANGLE 

V/VINF 

RHO/RHOINF 

T/TINF 

P/P INF 

. 1 

l.tV'1 

.( 349 

*9995 

.0117 

• 0003 

88.4170 

• 0103 

4.2264 

22.3311 

94.3802 

2 

1.C152 

• 1354 

2. Cl 46 

. 1179 

.3159 

48.3361 

.0245 

4.2230 

22.3206 

94.2596 

3 

1.030?' 

• C 360 

l.f-296 

.0190 

• C31 5 

32.C664 

• 0386 

4.2175 

22.3016 

94.0567 

4 

1.L452 

.( 365 

l.( 446 

.0194 

• 0486 

21.7493 

•0526 

4.2099 

22.2742 

93.7714 

5 

1.0603 

.( 37L 

1.1 597 

..•'134 

• C63 8 

16.1277 

.0667 

4.2001 

22.2384 

93.4034 

6 

1.0754 

.i 375 

1.C747 

.0101 

.0783 

12.9421 

.0608 

4.1302 

22.1941 

92.9527 

7 

2 ,C9t4 

.'*38- 

l.f 898 

.*179 

• 0934 

10.8570 

• 0946 

4.1740 

22.1415 

92.4194 

8 

1.1255 

.L3&6 

1.3048 

. >177 

.1079 

9.3428 

• 1066 

4.1577 

22.0806 

91.8036 

9 

1.1205 

.0391 

1.1199 

.0176 

• 1221 

8.1078 

.1228 

4.1390 

22.0114 

91.1056 

U 

1.1356 

.1 396 

*.1349 

, 3174 

• 1361 

7.2955 

•1368 

4.1181 

21.9339 

90.3258 

11 

1.1507 

• ( 4C 2 

1.1499 

.3173 

.1501 

6.5784 

• 1508 

4.0948 

21.8482 

89.4648 

22 

1.; 657 

41 7 

1.1650 

.0172 

.1639 

5.9916 

• 1646 

4.0692 

21.7542 

88.5233 

23 

2.1806 

.141 2 

1.1820 

. >171 

.1778 

5.4997 

.1787 

4.0413 

21.6521 

87.5019 

14 

1.1958 

• C 427 

-.1951 

• U7> 

• 1916 

5.0828 

.1926 

4.0109 

21.5419 

86.4017 

15 

1.2179 

• C 4? 3 ' 

1.2101 

.3170 

• 2054 

4.7228 

• 2065 

3.9780 

21.4235 

85.2237 

16 

1.2259 

.(428 

1.2252 

.2269 

.2193 

4.411J 

• 2204 

3.9428 

21.2970 

83.9689 

17 

1.2410 

• C43 3 . 

1.2402 

.0269 

• 2334 

4.1354 

• 2343 

3.9U50 

21.1625 

62.6387 

28 

2.256C 

42 6 

1.2553 

.0169 

• 2474 

3.9001 

• 2481 

3.6646 

21.0200 

81.2346 

19 

: .2711 

• l 444 

1.27(3 

• 0269 

.2621 

3.6835 

.2620 

3.6216 

20.8695 

79.7581 


ANGULAR 

LOCATION NO. 

3# AT THETA 

- 6.0031 

OEGREES 







I 

RP/D 

R/D 

X/D 

VR/VINF 

VX/VINF 

FLOW ANGLE 

V/VINF 

RHO/RHOINF 

T/TINF 

P/PINF 

1 

i.CCll 

.1(46 

.9957 

.0350 

.0029 

85.2522 

.0409 

4.2139 

22.2977 

93.9609 

2 


• )Uc 

l.Cx27 

.7543 

• 0185 

71.1386 

• 0525 

4.2102 

22.2745 

93.7792 

3 

1.0314 

• it:7e 

l. f . 257 

.“585 

.0367 

58.4397 

• 0642 

4.2045 

22.2454 

93.5316 

4 

1.0465 

• 1C94 

1.04C8 

.3573 

.5516 

48.C188 

.0760 

4.1970 

22.2102 

93.2169 

5 

2.0 616 

•me 

1 . 05 56 

.0553 

iC67C 

39,4164 

.0878 

4.1875 

22.1669 

92.0335 

6 

1.C767 

.1125 

1 .f 708 

.•?541 

• i8i 7 

33.5384 

.0997 

4.1761 

22.1212 

92.3799 

7 

1.0918 

.1161 

1.0639 

• 0536 

• C963 

29.0690 

• 1117 

4.1625 

22.0670 

91.0546 

8 

1. 7.069 

.1157 

i.lv ; 9 

• 3531 

.1107 

25.6199 

• 1239 

4.1469 

22.0058 

91.2559 

9 

3 .1221 

.1173 

1.1159 

• ibZb 

.1249 

22.8317 

.1362 

4.1291 

21.9376 

90.5023 

10 

1.1372 

.2164 

1.1309 

.0522 

.1389 

20.5983 

• 1487 

4.1)91 

21.8619 

09.0321 

11 

1.1522 

. 12( 4 

1.2460 

.>519 

• i529 

18.7485 

• 1613 

4.0868 

21.7784 

09.0035 

12 

1.1674 

.1221 

1.1610 

. >516 

.1667 

17.2C44 

.1741 

4.0622 

21.6866 

88.0951 

13 

1 . 1 62 5 

.1236 

1.1760 

.0514 

• 1806 

15.8878 

.1871 

4.0352 

21.5863 

07.1050 

14 

1.1976 

• 2 252 

2.1911 

• 5312 

.1944 

14.7565 

.2004 

4.0058 

21.4768 

06.0317 

15 

1.2128 

• 1 266 

1.2C61 

• 0510 

.2083 

23.7666 

• 2139 

3.9739 

21.3576 

04.0737 

26 

2.2279 

.1263 

1.2211 

.0509 

• 2222 

12.8966 

.2276 

3.9395 

21.2282 

03.6293 

17 

2.2 43r 

.1249 

a . 2362 

.35*7 

• 2364 

12.1097 

• 2416 

3.9025 

21.0801 

02.2972 

16 

1.2531 

325 

1.2612 

.55w6 

.2504 

11.4144 

.2559 

3.0629 

20.9364 

80.0760 

19 

1.2732 

.1331 

2.2662 

. 0 5i‘ 3 

.2652 

10.7386 

• 2705 

3.8206 

20.7726 

79.3645 
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APPENDIX A 


ANGULAR LOCATION ND. 23* AT THETA ■ 86.^.*.- DiGRc-S 


I 

RP/Q 

R/D 

x/o 

VR/VINF 

VX/VI NF 

FLO* ANGLE 

V/ViNF 

RHO/RN JINF 

t/tinf 

P/PINF 

1 

1.309 2 

l . 3061 

.T9-3 

.31*5 

.62*3 

26.4427 

.6975 

1.6544 

11.9337 

19.7001 

2 

1.3464 

1.3451 

• C- 9 4 1 

• 31 33 

.6233 

26.6613 

.6974 

1.7271 

11.9569 

20.6503 

3 

3 .3t)75 

.1.3*41 

.0468 

.3124 

.6248 

26.5626 

.6993 

1.8 >11 

11.9305 

21.4079 

A 

1 .4256 

1.4222 

995 

.3133 

.6251 

26.5616 

. 70CG 

1.6770 

11.8797 

22.2965 

5 

1.4658 

1.4(22 

.1022 

.3140 

• 6233 

26.554? 

.7026 

1.9554 

U.803t 

23.0794 

6 

1.5v49 

1 .5'.. 3 

• C '.54 

.2154 

.0313 

26.5492 

• 7059 

2.0367 

11.7052 

23.8401 

7 

3 .5441 

1.5‘i 3 

.1 . 77 

.3172 

• o348 

26.5527 

. 7C9 7 

2.12 15 

11.5070 

24.3017 

e 

1. 583 2 

1.5743 

.1134 

.3194 

.0386 

26. 5630 

.7142 

2.2112 

11.4502 

25.3076 

9 

1 .6223 

1.6384 

.1132 

.3219 

.0433 

26.5812 

.7193 

2.3035 

11.2962 

26.0208 

10 

1.6615 

1.6574 

.1159 

.3246. 

• 6431 

26.6f4r 

.7246 

2.4117 

11.127Q 

26.7242 

11 

;.?< 16 

1 .69t5 

» 1 1 8 j 

.3274 

.6532 

26.6338 

.7307 

2.5155 

10.9441 

27.4207 

12 

1.7297 

i.73‘ 5 

.12-4 

• 33%a 

• 6586 

26.6691 

.7369 

2.6153 

10.7494 

20.1133 

13 

1.7739 

1.7744 

.:;4i 

.3341 

• 6641 

26.7088 

.7434 

2.7317 

10.5446 

28.8040 

14 

2 ,SZ9C 

1.6236 

.2268 

.3376 

.6693 

26.7519 

.7501 

2.8551 

10.3317 

29.4906 

15 

1.6571 

1.8526 

.1295 

• 342 2 

.6755 

26.798C 

.7569 

2.9862 

10.1126 

30.1980 

14 

1.8953 

1.6917 

• 1323 

. 3449 

• 0812 

26.8469 

.7637 

3.1253 

9.889V 

30.9071 

17 

1.9254 

1.4?l 7 

• 1?50 

• 3485 

. 68 72 

26.8922 

.7706 

3.2730 

9.6641 

31.6306 

16 

.974 5 

3 .9697 

.13 77 

.3522 

.6931 

26.9339 

.7774 

3.4299 

9.4399 

32.3742 

19 

2.C137 

2.ltfcd 

.1405 

• 35tO 

.6937 

26.9992 

.7841 

3.5964 

9.2161 

33.1445 


ANGULAR 

LOCATION NO, 

,24, aT THETA 

• 9C.OOOO 

DEGREES 







I 

RP/Q 

R / D 

x/ z 

VR/VIN* 

VX/VINF 

FLO* ANGLE 

V/VINF 

RHQ/RHOINF 

T/TINF 

P/PINF 


1.3491 

1 .3491 


.2924 

.6553 

24.T135 

.7185 

1.5205 

11.3201 

17.2127 

2 

1.3922 

1.2922 

-. f- ■} 

• 29b6 

. 6543 

24.2941 

• 718C 

1.5937 

11.3353 

16.0653 

3 

1 .4352 

i.4352 

-.LC '‘"I* 

.29 7:* 

.6544 

24. 4 .*8 2 

• 7165 

1.6685 

11.3199 

16.0870 

4 

3 .4783 

1.4783 

-.C'.'Q 

.2987 

• 6 550 

24.5124 

.7199 

1.7456 

11.2760 

19.6836 

5 

1.5214 

1.5214 

-•OC ’ - 

.3G.-7 

• 6554 

24.6*19 5 

• 7222 

1.6255 

11.2059 

20.4563 

6 

1.564;* 

1 . £ 645 

-.o:oo 

• 30 30 

,6538 

24.6965 

.7252 

1.9.00 

11.1114 

21.2095 

7 

1 . 6(-7 6 

A.ti 76 

-•r.. 03 

• 3056 

.6618 

24.7844 

.7290 

1.9961 

10.9947 

21.9464 

5 

1.6576 

l.<5t6 

-.0' .. f . 

.*065 

• 6634 

24.8727 

• 7334 

2.0885 

10.8577 

22.6 704 

9 

1.6937 

1.6927 

-.ccoo 

.3117 

.6695 

24.9622 

• 7384 

2.185;* 

10.7025 

23.3846 

16 

:.?J66 

• .7368 

. «• 

.315* 

.6739 

25.0535 

.7438 

2.2877 

10.5309 

24.0920 

11 

1.77Q9 

i . 7799 

-.L 

• 3196 

.6737 

25.1466 

• 7496 

2.3968 

10.3451 

24.7955 

12 

1.623C 

1.623C* 

-•CC-QQ 

.3223 

.6837 

25.2414 

.7558 

2.5129 

10.1471 

25.4981 

13 

I .866' 

1.866'- 

— • L v L 0 

• 3262 

.6889 

25.3375 

.7622 

2.6364 

9.9391 

26.2029 

14 

1.9091 

1.9L91 

. -.eocc 

.3352 

.6943 

25.4343 

• 7608 

2.7679 

9.7233 

26.9135 

.15 

1 .9522 

1.9522 

-.coco _ 

. 3342 . 

.6997 

25.5325 

.7755 

2.9^02 

9.5020 

27.6336 

'.16 

1.9952 

1.9953 


.3383 

.7052 

25.6292 

• 7823 

3.0577 

9.2776 

20.3661 

1? 

2. 0384 

2.02*4 

.; >*;*> 

. 3424 

.71 v 7 

25.7247 

.7090 

3.2170 

9.0525 

29.1223 

19 

2 .v 614 

2.1 81 4 


.3435 

.7162 

25.8163 

.7956 

3.3868 

8.8294 

29.9034 

19 

2.1245 

2.12*5 

• «. • « 1 

.35- 3 

.7216 

25.9.-33 

.6020 

3.5o76 

8.6138 

30.7195 


FL OR FIELD VALUES FROI MARCHING CALCULATION 


ADDITION* 

L AXIAL 

LOCATION NC. 

1» AT X/D - 

-.1302 





I 

R/0 

VR/VINF 

VX/VINF 

FLOW ANGLE 

V/VINF 

RHQ/RHOINF 

T/TINF 

P/PJNF 

• i 

1.3623 

• 2870 

• 6670 

23.28C2 

.7261 

1.4791 

11.0865 

16.3985 

2* 

:.4i,6C 

.24*6 

. 66 j 3 

23.6535 

.7242 

1.5463 

11.1463 

17.2359 

' 3 

1.4497 

.242 3 

.6635 

22.7733 

.725 v 

1.6289 

11.1204 . 

18.1135 

4 

1 .*933 

• 29*1 

.6639 

23.3901 

.7261 

1.7062 

11.0855 

18.9140 

5 

1.537C 

.24t5 

• 6o*»7 

2**. . 394 

.7272 

1.7885 

11.0321 

19.7305 

6 

1.5807 

.2992 

.66ofc 

24.1753 

.7306 

1.6736 

I0.945r 

20.5169 

7 

2 .624? 

.•I-.22 

.6692 

24,3021 

.7343 

1.9621 

10.8313 

21.2519 

0 

2 ,ttiK 

• 3i L'4 

.6 725 

24,4241 

. 7326 

2.0550 

1C. 6965 

21.9810 

9 

1.711? 

• 3( ot* 

• 6762 

24.5433 

.7434 

2.1531 

LI. 5439 

22.7019 

1 J 

1.7552 

• 2*i< 4 

• 6?0 4 

24.6629 

.7437 

2.2571 

10.3757 

23.4190 

1* 

j .799.* 

.3:67 

. 6c 49 

24.781? 

.7543 

2.3677 

1C .1939 

24.1363 

12 

1.8427 

• 320 1 

.6896 

24.90*2 

• 7603 

2.4856 

10.0036 

24.8574 

13 

1,0363 

.3242 

.6946 

25.1183 

.7685 

2.6113 

9.7979 

23.5954 

14 

2.93:*C 

.3203 

.6948 

25. *355 

.7730 

2.7454 

9.5877 

26.3217 

15 

1.9727 

.332 5 

• 7C5C 

c5. 21*35 

.7795 

2.8882 

9.3715 

27.0671 

1 6 

2.**:73 

. 3 2 6 n 

.7: 3 

26.2662 

.7861 

3.0407 

9.1511 

27.8254 

17 

2 • C 610 

.34. , 

.7:56 

25.4784 

.7927 

3. 2035 

6.9233 

28.6319 

10 

2.1:47 

.1452 

.7 2?,9 

25.587') 

.7993 

3.3775 

0.7051 

29.4012 

19 

2.1403 

.2442 

• 72 61 

? 5 . b8 29 

.8)57 

3.5617 

0.4847 

30.2195 


ADDITIONAL AXIAL 

L0CAT1CN NO. 

2# A r X/D - 

-. .674 





I 

F/0 

VR/VirF 

VX/VINF 

FLOW ANGLE 

V/VINF 

RHQ/RHOINF 

T/TINF 

P/P IMF 

1 

1.37->2 

.2011 

.6773 

22.5497 

• 733i 

1.4358 

10.8690 

15.6058 

2 

1.4194 

• c o4tl 

.6728 

22.9425 

.73 06 

1.5056 

10.9463 

16.48D4 

3 

2.4636 

.2872 

.67-7 

23.1539 

.7306 

1.5925 

10.9472 

17.4339 

4 

1.5C79 

.2^96 

.6717 

23.3241 

.7314 

1.6709 

10.9198 

18.2457 

5 

1.552? 

.29 cc 

. 6? 24 

23.49L2 

. 7332 

1.7523 

10.8650 

19.0366 

A. 

1.5964 

. ?9t4 

.6 73 8 

23.6725 

• 73?7 

1.6386 

10.7852 

19.8318 

7 

1.6407 

.2467 

• 6?0l 

23.3391 

.739 2 

1.9287 

10.6780 

20.5943 

8 

1.6P5C 

• ;L2? 

.6790 

23.9928 

.7433 

2.r229 

10.5402 

21.3378 

9 

1.7292 

.3059 

• to 25 

24.1426 

.7479 

2.1224 

1C. 3990 

22.0720 

ID 

2.^735 

• 3' 98 

.6864 

24.2380 

.7531 

2.2278 

10.2352 

22.8024 

U 

:.8:78 

.2236 

.69 .7 

24.4319 

.7566 

2.3401 

10.0566 

23.5334 

32 

1.062C 

.2179 

.6952 

24.5717 

.7645 

2.4598 

9.8661 

24.2687 

13 

1.9L63 

.2221 

. 7 U C 5 

24.7105 

.7700 

2.5875 

9.6661 

25.0112 

14 

1.9536 

.22/4 

.7- 44 

24.3474 

.7759 

2.7238 

9.4583 

25.7627 

15 

1.994e 

. 33C6 

.72 uO 

24.9817 

.7833 

2.8692 

9.2445 

26.5241 

16 

2.3391 

.2252 

.7151 

25.1125 

.7898 

3.0244 

9.0263 

27.2995 

17 

2.;. *54 

.3395 

.721 3 

2J.239'* 

.7953 

3.1904 

8.0059 

23.C944 

18 

2.1276 

.2438 

.7254 

25. 3608 

• 802 3 

3.3678 

8.5853 

28.9135 

10 

j* . > T ■ O 

.UK 

.71 :S 

✓4.44QA 

. AA9 A 

3. 5855 

8.3663 

29.7463 
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APPENDIX A 


ADDITIONAL 

, AXIAL 

L GC ATI OK NO. 46/ 

AT X/0 - 

-4.6689 





X 

R/0 

VR/VINF 

vx/vinf 

FLOW ANGLE ~ 

V/VINF 

.RHO/RHOINF 

T/TINF 

P/PJNF 

1 

2.7691 

• Cc73 

.4043 

• <60 7 

.9044 

• 4326 

4.8848 

2.1132 

2 

1.9171 

•i 13"5 

• 9:7a 

.3494 

.9077 

.4504 

4.7561 

2.1421 

3 

2.C649 

•L?r c 

.9212 

1.2599 

• 9114 

• 4728 

4.6136 

2.1813 

A' 

2.2128 

•C272 

.9144 

1.7062 

.9148 

• 4989 

4.4810 

2.2358 

5 

2. stl, b 

•i 354 

.9175 

2.2C76 

.9182 

.5317 

4.3473 

2.3116 

6 

2.5035 

.C447 

.9203 

2. 780o 

.9214 

• 5726 

4.2207 

2.4167 

7 

2.6554 

•r 563 

.9228 

3.4322 

.9244 

• 6236 

4.1024 

2.5584 

8 

2.8o43 

• t 572 

.9248 

4.1552 

.9272 

• 6871 

3.9912 

2.7422 

9 

2.9522 

.C6C0 

.9255 

4.933r, 

.9299 

.7652 

3.8848 

2.9725 

10 

3. 10SC 

• L 934 

.9279 

5.7464 

.9325 

.6603 

3.7812 

3.2531 

11 

2.2479 

• 1C 71 

.9289 

6.3772 

.9351 

.9753 

3.6795 

3.5387 

12 

3.3958 

. 1 2C9 

. 4297 

7.4115 

.9376 

1.1128 

3.5812 

3.9350 

13 

3.5437 

.1347 

.9302 

8.2382 

.9399 

1.2753 

3.4685 

4.4490 

14 

3.t9! 6 

• I4e2 

.9312 

9.0545 

.94*9 

1.4658 

3.4052 

4.99X5 

15 

3.0394 

.1614 

.9299 

9.3443 

.9438 

1.6855 

3.3316 

5.6155 

16 

2.9872 

.1743 

.9291 

1..-.5241 - 

• 9452 

1.9374 

3.2740 

6.3429 

17 

4.1352 

.26# 2 

.9279 

•1.3466 

• 9464 

2.2159 

3.2265 

7.1496 

18 

4.2831 

• 1665 

• 9256 

12.1035 

.9468 

2.53C9 

3.2082 

8.1195 

19 

4 . 431v’ 

• 2<it9 

• 9 24U 

12,7372 

.9473 

2.8557 

3.1897 

9.1009 


ADDITIONAL 

AXIAL 

LOCATION NO. 46/ 

at tn • 

-4.9275 





I 

R/0 

VR/Vlt*F 

VX/VInP 

FL3« ANGLE 

V/VInF 

RH0/RH01HF 

T/TINF 

P /P INF 

1 

1.771C 

• l • r 3 

• 9..5S 

• 40J 6 

.9055 

.4269 

4.8417 

2 • 0o69 

2 

1.9237 

• C 12a 

.9009 

.7619 

.909*4 

• 4444 

4.7C7C 

2.1.917 

3 

2.v 765 

• x 160 

.9126 

1.1776 

• 9128 

.4663 

4.5 58 3 

2.1256 

4 

2.2293 

•i 2io 

«916v 

1.6C32 

• 9134 

.4918 

4.4183 

2.1731 

5 

2.362C 

• C 334 

.9194 

2.0776 

.420 0 

.5238 

4.2764 

2.2398 

6 

2.5348 

••‘422 

.9225 

2.5294 

• 9234 

• 5635 

4.1412 

2.3334 

7 

2.687f 

.Cal 4 

.9252 

3.2396 

.9256 

.6131 

4.C151 

2.4al6 

8 

2.3403 

• L62o 

.9274 

3.9352 

• 4246 

.6752 

3.8974 

2.0314 

9 

2.9931 

• i 7t 3 

• S29 j 

4.6915 

• 9324 

.7519 

j 3.7663 

2.8468 

10 

3.1456 

• C094 

.9303 

5.4892 

.9351 

• 6458 

3.6798 

3.1124 

11 

2.293? 

.1030 

.9320 

c. 3093 

.9377 

.9596 

3.577C 

3.4325 

12 

3.4512 

.116 6 

• 93 28 

7.1360 

.9401 

1.C960 

3.479C 

3.8131 

13 

2.604C 

.13C5 

.9333 

7,9541 

.9424 

1.2576 

3.3679 

4.2006 

14 

3 . 7 5fc e 

.J4AU 

.9333 

9.7735 

.9444 

1.4472 

3.3073 

4.7862 

15 

3. 9095 

• * 572 

.933C 

9.5633 

• 94o* 

1.6659 

3.2371 

5.3925 

lb 

4.0623 

.1701 

.4320 

1C. 3447 

• 9474 

1.9168 

3.1836 

6.1,23 

17 

4.2151 

.2 o2i 

.9309 

11.0702 

.9435 

2.1937 

3.1406 

6.8896 

18 

4 . 367 d 

.1945 

.9297 

11.9200 

.9488 

2.5C66 

3.1272 

7.6393 

19 

4.5216 

.2i 49 

.9268 

12.4041 

.9492 

2.8289 

3.1137 

8. 0981 


ADDITIONAL 

AXIAL 

LOCATION NO. 

59/ 

AT A/0 « 

-9.0204 




1 

I 

P/0 

VR/VINF 


VX/VINF 

FLOW ANGLE 

V/VINF 

RHQ/RHOINF 

T/TINF 

P/PINF 

1 

1.7848 

• Cl C 9 


.9196 

• .577 

.9196 

.3564 

4.2931 

1.5302 

2 

2.L22C 

• CCc3 


.9245 

.3892 

.9245 

.3738 

4.0994 

1.5322 

3 

2.2591 

• om 


.429o 

.6816 

.9297 

.3953 

3.8951 

1.5396 

4 

2.4953 

• ( 117 


.9348 

.96 35 

.9349 

.4208 

3.6854 

1.5508 

5 

2.7235 

•f. 2f 3 


.94.2 

1.2367 

.9434 

.4519 

3.4675 i 

1.5671 

6 

2.971.6 

• X.25C 


.9454 

1.5163 

.9453 

.4890 

3.2515 ! 

1.5912 

7 

5.2» 78 

• i'3l u 


.9:1 5 

1.3183 

• 9529 

• 53 34 

3.0425 ' 

2.6229 

8 

3.4450 

.1215 


.9tax 

2.1318 

.9557 

.5667 

2.8474 

1.6705 

9 

3.6821 

.0421 


• 9590 

2.5163 

• 9600 

• 6516 

2.6742 \ 

1.742b 

1C 

2.9193 

•L 5L 2 


.9 o22 

2.9890 • 

.9635 

.7321 

2.5288 

1.6514 

11 

4.1555 

• C‘*S6 


• 9c 45 

3.5531 

.9663 

.8318 

2.4129 

2. 1.-69 

12 

A. 3936 

.i K 7 


.9639 

4,1939 

.9635 

.9536 

2.3246 

2.2266 

13 

4. 6308 

• 1.823 


.9655 

4.3797 

.9711 

1.C988 

2.2590 

2.4031 

14 

4 . 8 66 C 

.0945 


.9654 

a • 6L ; %. 2 

.971, 

1.2693 

2.218C 

2.0154 

15 

5.1:51 

• li-14 


.9658 

6.3153 

.97:6 

1.4622 

2.1942 

i 3.2084 

16 

5. 3423 

.1166 


.9643 

,7.1594 

.97*6 

1.6837 

2.1949 

\ 3.689C 

17 

? . 5 79 ? 

• 13C2 


.9625 

7.75x3 

,9713 

1.9134 

2 . 206b 

4.2225 

18 

5.6106 

• i42 4 


• 9590 

8.4982 

.9703 

2.1720 

2.2468 

4.60CO 

19 

6.C53S 

.1527 


.9573 

9.1393 

.4694 

2.4304 

2.2863 

5.5565 

ADDITIONAL 

. AXIAL 

LOCATION NO, 

>60/ 

AT X/0 - 

-1C. >903 




\ 

I 

R/D 

VR/VINF 


VX/VINF 

FLOW ANGLE 

V/VINF 

RHO/RHDINF 

T/TINF 

P/PINF 


. * 7854 

• oCCa 


.9209 

.0476 

.4209 

.3501 

4.2424 

1.4654 

2 

2. * , .329 

•( V LI 


• 92 59 

• 3733 

.9254 

• 3676 

4.0439 

1.4966 

3 

2.2653 

• c;* 7 


.9511 

.6598 

.9312 

.3894 

3.8351* 

1.4935 

4 

2.527? 

•C i i 2 


.9365 

.9 330 

.9366 

.4154 

3.62C3 

1.5039 

5 

2.7752 

•f 2 4 7 


.9419 

1.1985 

.9421 

.4472 

3.3972 

1.5191 

6 

3.C227 

.(242 


.9473 

1.4646 

.9476 

.4351 

3.1763 

1.5427 

7 

3.27;<? 

.< 2t9 


.9:24 

1.7406 

.9529 

.5302 

2 «9o28 

1.5710 

8 

3.5175 

•1341 


.4372 

2. •>415 

.9578 

.5841 

2.7637 

1.6143 

9 

3.765C 

• C4C2 


.4613 

2.2948 

.9621 

.6491 

2.5869 

1.6793 

10 

4.C124 

•L 47? 


.4645 

2.8310 

.4657 

.7289 

2.4391 

1.7777 

11 

4.2599 

• i 5e 6 


.9360 

2.3631 

.9635 

.8273 

2.3227 

1.9215 

12 

4.5C73 

.f 672 


.9103 

3.9775 

.9716 

,9475 

2.2360 

2.1186 

13 

4.7546 

• C7t5 


.9689 

4.e429 

.9721 

1*0909 

2.1742 

2.3726 

14 

5.? » 22 

• ( 4i 4 


• 9668 

5.3477 

• 9730 

1.2587 

2.1365 

2.689* 

15 

L. 2496 

.1022 


• 9o31 

d.0493 

.9735 

1.4470 

2.1171 

3.2653 

16 

5.4971 

• 1**4 


• 9666 

c. 7019 

.9733 

1.6615 

2.1223 

3.5262 

17 

5.744? 

.1257 


• 9c46 

7.4643 

.9730 

1.8681 

2.1365 

4.0377 

16 

6.9920 

.1377 


.9621 

a. 20*3 

.9719 

2.1395 

2.1816 

4.6679 

19 

6.2394 

. 24fcc 


.9595 

6.0371 

.9709 

2.3901 

2.2247 

5.3173 


Figure A.2(£) - continued 
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APPENDIX A 


STRsANLINc TRA JECT3RT CALCULATION 


A3 STREAMLINES CALCULATED 


STREAMLINE NO. 1# STARTING AT X/O • 1.2753# A/O - ,1:444 

(CORRESPONDS TO T4ETA • 2.0C DEGREES# AP/O - 1.2711) 


X/D 

A/D 

1*2703 

. C 4 4 4 

1.2413 

. 4 6 £ 

1.2153 

.C 433 

a.: 9'-' 3 

• ibl £ 

1.1653 

.0552 

1.1403 

• 06DC 

1.115; 

•>.c7v 

1.0903 

. 373C 

1.5653 

.-■933 

l.*5« 3 

.i.232 

1.0153 

.1668 

.4 9:’? 

.263a 

.96 5 3 

.3498 

.9353 

.4353 

.91^ 1 

.496*. 

.885 3 

.3 5‘ 9 

• 3 60 3 

• 6 cu 2 

.3353 

.645*. 

.0103 

.66o9 

.7853 

.7257 

.7 6 -.' 3 

.762 

.7353 

.7961 

.71 3 

• 3264 

.585 3 

. .6539 

.6603 

.5679 

.63(3 

.920 9 

.605 3 

.9471 

.5 0-: 3 

.9721 

.5553 

.9961 

.5303 

A •**" 1 9 a 

.505 3 

1.0412 

.4803 

... 624 

.4553 

1.0629 

.4 3. 3 

1.1C26 

• 4( 53 

2 • .216 

• 3 6 C 3 

1.140u 

.2553 

...577 

.3253 

1.1762 

.3(03 

1.1947 

.2753 

2.2106 

.2503 

1.2261. 

.2253 

1.24C9 

• 2i<>3 

..1554 

.1753 

1.26*4 

• 1503 

1.263A 

.1253 

. 1.2962 

.1003 

1.309C 

..753 

2.3214 

.3 5:2 

2 • jS34 

.3253 

1. 3451 

-.;r*7 

2.3530 

-.2161 

i.*.422 

-.4276 

1.5C77 

-.0003 

..5S3‘: 

-.6202 

1.5939 

-1 .vS* t 

2 .oofli 

-1.2129 

1. o543 

-1.4243 

1.6773 

-1.5.' 5 5 

.« o93c 

-1.3 : 70 

1.7C93 

-2 — 22 4 

1.7223 

—2.2(96 

1.7316 

-2.421C 

1. 7407 

-2.6325 

1.74 8 4 

-2.3137 

1.754. 

-3.0251 

1. /5 97 

-3.2( 64 

1.7639 

-3.4178 

1.7681 

-3.6292 

i.77ia 

—3 .61 ^ 4 

1.7740 

-4,0219 

1.7774 

-4.2333 

1.7799 

-4.4195 

1.7619 

-4.6259 

1.7839 

-4.3f 72 

1.7854 

-4.9175 

1. 7063 


Pigure A. 2(g) - continued 
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APPENDIX A 


STREAMLINE NQ.42, STARTING AT X/D ■ -4.3769, R/o ■ 4.3233 

X/D K/Q 

-4.3769 4.3233 

-4.5663 4.3717 

-4.7996 4. ‘,187 

-4.9176 4.4442 


STREAMLINE NO. 43, STARTING AT X/O • -4.6915, R/O - .4.4391 

X/D R/O 

-4.6915- 4.4391 

-4.9029 ' 4.4861 

-4.9175 4.4693 


MAGNETIC FIElO COMPONENTS 


ANGULAR 

LOCATION 

NO. 

1, AT THETA 

- t.:*c*:o 

DEGREES 






I 

RP/D . 


8/B1KF. 

3 / 6 1 NF 

B/9INF 

• 8-ANGLE 

B/BINF 

0X/BINF 

8Y/BINF 

BZ/BINF. 




(PARALLEL ) 

(PER?) 

(IN-PLANE) 

- (IN-PLANE) ^ 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT 

1 

1 .0030 


c.tccc 



> 





2 

i *115.'’ 


.0693 

19.5068 

13.7949 

r -89.7964 

4.2200 

.0347 

-9.7544 

2.984C 

3 

1.6 31>1 


• 1292 

15.7355 

11.63 21 

- -84.5076 

4.2156 

.0646 

-7.5177 

2.9809 

4 

1 . 46 1 


.1898 

11.7087 

8. 2804 

• -39.0712 • 

4.2089 

• 9949 

-5.0544 

2.976* 

5 

6.. 2 


.25i 9 

9.6738 

t. 9841 

• -88.5445 

4. 1990 

.1254 

-4.9369 

2.9697 

6 

1.0752 


* 3121 

3.5656 • 

6. .*615 

-87.9137 - 

4.1383 

.1560 

-4.2833 

2.9616 

7 

1.0903 


.3731 

7.5d05 

5.3667 

-87.1021 • 

4.1744 

• 1866 

-3.7903 

2.9518 

8 

I.IC‘3 


.4337 

6.843* - 

4^3521 

-96.3759 

4.1582 

.2169 

-3.4241 

2.9403 

9 

1.1204 


• 4436 

5 • 2 9 ; 3 

4.4651 

• -05.5*oo 

4.1395 

.2466 

-3.1476 

2.9271 

10 

i • 1 35 4 


.5526 

1.6932 

4.id4* 

-34.6424 

v.:i85 

.2762 

-2.9456 

2.9122 

n 

1.1 >•’* 6 


•el(l 

5.5159 

3.9241 

- -83.6888 

4.0951 

.305? 

-2.7589 

2.8957 

12 

7—1655 


• 6C6l '■ 

5.2i 39 

2.7167 

-82.7202 

4.C69 3 

• 3330 

-2.6169 

2.8775 

13 

1.1655 


.721 1 

4.9412 

3.3301 

-81.7056 

4.0412 

.3601 

-2.473k 

2.8575 

14 

1 • 5 95 1 


. 7 722 

4.7254 

2.3879 

-80.7244 

4.0106 

• 3861 

-2.3p42 

2.8360 

15 

1.2106 


.6219 

4.5212 

2.2444 

-79.6968 

3.9777 

• 411C 

-2.2606 

2.8127 

16 

1.2257 


• c 64L 

4.3214 

3. *161 

-78.6277 

3.9424 

.4345 

—2**632 

2.7877 

17 

} .24:>7 


.9132 

4.J o99 

3.0185 

-77.6478 

3.9048 

.4566 

• -2.9850 

2.7611 

16 

1.2 5a 8 


.9543 

4.*1 79 

2.4202 ' 

• -76.0392 

3.8647 

.4772 

-2.9*91 

- 2.7320 

19 

1.2 7'.- 6 


1 .v Lt l • 

3.82.9 

2.7928 

-75.3336 

3.8209 

• 5000 

-1.9*04 

2.7018 

ANGULAR 

LOCATION 

NO. 

2# AT THETA 

• 2.::u 

DEGREES 






I 

*P/D. 


3/6INF 

3/BINF 

8/3INF 

' 8-ANGLE 

B/BINF 

3X/BINF 

8Y/8INF 

BZ/BINF 




( PARALLEL) 

(PERP) 

( IN-PLANE) 

• UN-PLANE) 

(NORMAL) 

(RESULTANT) 

(RESULTANT ) 

(RESULTANT 

1 

l.OCM 


• C 437 








2 

i.C 152 


.10 3 3 

19,6130 

13.3129 ' 

"-91.2200 

6.8133 

-.208 C 

-9.7650 

4.8177 

3 • 

:.03f2 


* I 626 

I4.V7J3 

1L. 525: 

-90.7384 

6.6524 

-.0959 

-7,4417 

4.7039 

4 

1.0453 


• 2 226 

1 1 . 69C * 

8.2066 

-90.0175 

6.4249 

-.cua 

-5.8131 

4.5431 

5 

l.l 6.-3 


• 2 faL 2 ' 

9,6647 

0.9237.- 

-89.2630 

6.2297 

• 0630 

-4.8954 

4.4051 

6 

;.r75 4 


.336 3 

9.5682 

6.177: 

* -88.4134 

6.0036 

.1172 

-4 .2 *60 

4.2452 

7 

1.0904 


.39*8 

7.*a40 

5.3160 

-87.5814 ' 

5.7295 

.1587 

-3.7571 

4.0514 

6 

i . 1 06 : 


.4621 

6 • 05 iC 

4.3053 

-8o.fc 891 

5.5336 

• 1962 

-3.3922 

3.9129 

9 

1.1206 


.508 4 

6. 3w59 

4.4193 

-35.7462 ■ 

5.3252 

• 2318 

—3.1167 

3.7655 

1C 

1. IS.'O 


.5635 

5.8963 

4.I35C 

' -84.8951 

5.1446 

• 2602 

-2.9123 

3.6378 

11 

1 . 15:7 


• 6175 

8.52*1 

3.8743 

' -34.C02C ' 

4.9735 

• 28o3 

-2.7246 

3.5168 

12 

:.i65? 


• t7i 4 

3 . 2 1 6w 

2.6645 

-83.1898 

4.6030 

.3073 

-2.5729 

3.3962 

13 

1 . 1 30 6 


.7221 

4.9419 

3.4755 

-02.3451 

4.6551 

.3276 

-2.4357 

3.2916 

14 

1.1958 


.7725 

4.7240 

3.3289 

-81.4203 

4. 5 r >30 

.3512 

-2.3275 

3.1841 

15 

: .21*9 


.82.5 

4.5131 

3.1872 

-80.4476 

4.3600 

.3740 

-2.2225 

3.0830 

16 

1.2239 


.6690 

4.3112 

3.0529 

-79.4017 

4.2217 

.3971 

-2.1219 

2.9852 

17 

1 .24U 


.9146 

4.2574 

2.9534 

-78.3632 

4*0025 

• 4212 

-2.1454 

2.6868 

16 

1 .2561 


• 9564 

4 • 0 j 4 9 

2.3555 

-77.28 J3 

3.9556 

• 4446 

-1.9696 

2.7971 

19 

1.2711 


1 .0011 

3. 8 2 v8 

2.7388 

-76.-275 

3.8218 

.4673 

—1 .8780 

2.7024 


Figure A. 2(h) - continued 
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ANGULAR LOCATION NO. 3» AT THETA • 6.0000 DE6REES 


l 

R P/0 

8/8XNF 

8/BINF 

8/8INF 

B-ANGLE 

B/BINF 

BX/BINF 

BY/B1NF 

3Z/BINF 



(PARALLEL l 

CPERP) 

(IN-PLANE) 

(IN-PLANE) 

(NQRNAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

l 

1.0011 

.1723 








2 

1.0163 

.2211 

19.7356 

13.8038 

-96.0993 

12.0951 

-.6977 

-9.7358 

8.5525 

3 

1.0316 

• 27t'b 

16.8020 

10.2979 

-93.6057 

1-1. 6519 

-.6326 

-7.2588 

8.2391 

6 

1.0665 

• 3188 

11.6097 

8.0369 

-92.16A3 

10.9525 

-.2126 

-3.6776 

7.7666 

5 

1.0616 

.3676 

9.8226 

6.7739 

-91.0505 

10.3597 

-.0876 

-6.7891 

7.3256 

6 

1.0767 

•6166. 

8.5379 

5.8698 

-89.9227 

6.6352 

.0036 

-6.1506 

6.1060 

7 

1.0918 

• 6651 

7.5758 

5.1965 

-88.9636 

7.5021 

• 0666 

-3.6725 

5.3068 

e 

1.1069 

.5138 

6.8516 

6.6877 

-87.9932 v - 

6.8321 

.1161 

-3.3127 

6.8310 

9 

1.1221 

.5626 

6.3036 

6.3063 

-87.0325 

6.3063 

• 1576 

-3.1639 

6.6592 

10 

1.1372 

• 61C8 

5.8859 

6.0167 

-86.1862 

5.9080 

• 1690 

-2.6339 

6.1776 

11 

1.1523 

.6592 

5.5098 

3.7568 

-35.2620 

5.5322 

.2196 

-2.6676 

3.9119 

12 

1.1676 

.7073 

5.2036 

3.5660 

-86.6222 

5.1768 

• 2637 

-2.6955 

3.6591 

13 

1.1825 

.7551 

6.9275 

3.3572 

-33.6896 

6.9066 

• 2692 

-2.3586 

3.6682 

16 

1.1976 

.8027 

6.7026 

3.2061 

-82.6689 

6.6555 

• 2899 

-2.2670 

3.2919 

15 

1.2128 

• 8699 

6.6867 

3.0578 

-81.7628 

6.6681 

• 3098 

-2.1399 

3.1653 

16 

1.2279 

• 8966 

6.2863 

2.9220 

-80.7752 

6.2662 

• 3312 

-2.0395 

3.0166 

17 

1.2630 

.9629 

6.3259 

2.8199 

-79.6669 

6.0980 

• 3563 

-1.9615 

2.8977 

18 

1.2581 

.9885 

3.9738 

2.7231 

-78.6669 

3.9595 

• 3850 

-1.8866 

2.7998 

19 

1.2732 

1.C336 

3.8123 

2.6216 

-77.1723 

3.8206 

.6116 

-1.8075 

2.7016 


ANGULAR LOCATION NO. 23# AT THETA • 86. COCO DEGREES 


I 

RP/D 

3 /B INF 

8/BINF 

B/3INF 

8-AN6LE 

8/BINF 

8X/BINF 

ar/aiNF 

BZ/BINF 



(PARALLEL) 

(PERP) 

( IN-PL ANE ) 

(IN-PLANE) 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

1 

1.3092 

1.1539 








2 

1.3686 

1.2165 

8.5172 

5.1726 

-169.3591 

11.6723 

-3.1669 

-1.8661 

8.3950 

3 

1.3875 

i.2577 

5.7915 

3.2113 

-165.2986 

7.1832 

-1.8891 

-1.2399 

3.0793 

6 

1.6266 

1 . 316c 

6.5566 

2.3G29 

-162.7691 

5.6296 

-1.2965 

-.9952 

3.9806 

5 

1.6656 

1.2738 

3.6260 

1.7525 

-138.9333 

6.8696 

-.9363 

-.6161 

3.6293 

6 

1.5069 

1.637 t 

3.3967 

1.6206 

-133.7266 

6.3286 

-.6962 

-.7258 

3. 0668 

7 

1.5661 

1.5(57 

3.K-93 

1.19C3 

-127.5731 

6.C165 

•>.5132 

-.6671 

2.8387 

8 

1.5832 

1.5787 

2.9117 

1.0300 

-120.5963 

3.7970 

-.3707 

-.6269 

2.6869 

9 

1.6223 

1.65o9 

2.7665 

•9114. 

-112.9630 

3.6526 

-.2513 

-.5931 

2.5828 

10 

1.6615 

1.7607\ 

2.6536 

• 6276 

-1W6.3051 

3.5698 

-.1666 

-.5669 

2.5101 

11 

1.7006 

1 .6307 \ 

2.5789 

.7682 

-95.5531 

3.685C 

-.0526 

-.5607 

2.6663 

12 

1.7397 

1.9272 

2. 52C2 

.7392 

-85.0703 

3.6306 

.0376 

-.5213 

2.6316 

13 

1.7789 

2.C3C7 \ 
2.1615 

2.6826 

.7373 

-76.28d9 

3.6192 

• 1236 

-.5-65 

2.6177 

16 

1.818C 

2.6593 

.7696 

-67.5293 

3.6158 

•2025 

-.6097 

2.6156 

15 

1.8571 

2.26C1 \ 

2.6526 

• 7S 76 

-59.3933 

3.6262 

.2836 

-.6793 

2.6227 

16 

1.8963 

2.3869 \ 

2.6656 

.8376 

-51.7332 

3.6525 

• 3666 

-.665C 

2.6613 

17 

1.9356 

2.5223 \ 

2.6666 

.0906 

-66.9236 

3.6900 

.6699 

-.6687 

2.6678 

18 

1.9765 

2.6666 

l 2.6689 

.9696 

-38.9013 

3.5389 

.5335 

— • 63C5 

2.5026 

19 

2.0137 

2.62C1 

' 2.6652 

1.0602 

-33.0877 

3.5966 

• 6162 

-.6015 

2.5630 

ANGULAR 

LOCATION 

NO. 26. AT THETA 

- 90 .0000 

DEGREES 






I 

RP/D 

B/BINF 

3/8INF 

B/BINF 

8-ANGLc 

8/BINF 

BX/BINF 

8Y/BINF 

BZ/BINF 



(PARALLEL) 

(PERP) 

(IN-PLANc) 

(IN-PLANE) 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

1 

1.3691 

1 .t 923 








2 

1.3922 

1.2663 

7.7661 

6.6839 

-151.7w32 

1C. 8080 

-2.9163 

-1.5700 

7.6626 

3 

1.6352 

1.1969 

5.3267 

2.9263 

-168.1508 

6.5661 

-1.7565 

— 1..912 

6.6630 

6 

1.6763 

1.2567 

6.1933 

2.0885 

-166.2260 

5.1963 

-1.1981 

-.0636 

3.6766 

5 

1.5216 

1.3166 

3.5221 

1.580? 

-139.7697 

6.6765 

-.8 53C 

-.7216 

3.1639 

6 

1.5665 

1.3866 

3.1363 

1.2765 

-133.7263 

6.0266 

-.6239 

-.6523 

2.6672 

7 

1.6076 

1.6532 

2.6766 

1 .3666 

-12 7. C 528 

3.7522 

-.6566 

-.6019 

2.0332 

8 

1.6506 

1.5316 

2.7CT1 

• 9232 

-119.2656 

3.5687 

-.3191 

— .5o95 

2.5236 

9 

1.6937 

1.6133 

2.5728 

• 8176 

-110.9357 

3.6686 

-.2666 

-•56J0 

2.6386 

10 

1.7368 

1.7C16 

2.6808 

.7696 

-lbl.2502 

3.37C7 

-.1036 

-.5197 

2.3636 

11 

1.7799 

1.7968 

2.6151 

.7039 

-91.6668 

3.3186 

-.0165 

-.6975 

2.3666 

12 

1.823C 

1.8992 

2.3690 

.6986 

-80.6905 

3.2956 

.0799 

-.6873 

2.3302 

13 

1.8660 

2.U95 

2.3626 

.6981 

-71.5558 

3.2926 

• 1562 

-.6682 

2.3261 . 

16 

1.9091 

2.1 281 

2.3295 

.7355 

-62.18/6 

3.3035 

•2627 

—.6605 

2.3359 

15 

1.9522 

2.2556 

2.3311 

.7869 

-56.2623 

3.3327 

• 3251 

-.6515 

2.3565 

16 

1.9953 

2.3920 

2.3315 

• 8636 

-66.9660 

3.3762 

.6G73 

-.6339 

2.3859 

17 

2.C386 

2.5383 

2.3372 

• 9186 

-60.6107 

3.6262 

• 6966 

-.6211 

2.6261 

18 

2.0816 

2.6966 

2.3669 

• 9966 

-36.6920 

3.6919 

.3793 

-.6010 

2.6691 

19 

2.1265 

2.8613 

2.3530 

1.0732 

-29.3216 

3.5676 

• 6ol7 

-.3716 

2.5226 


Figure A. 2(1) - continued 



APPENDIX a 


ADDITIONAL AXIAL 

LOCATION NO. 1. 

AT X/0 • 

-.0302 






1 

R/0 

8/BINF 

B/BINF 

B/BINF 

B-ANGlE 

B/BINF 

3X/9INF 

8Y/BXNF 

BZ/BXMF 

1 

1*3623 

(PARALLEL) 
1 . t 74^.- 

( PERP) 

( IN-PLANE) 

(IN-PLANE) 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

2 

1*4060 

1.1148 

7.5036 

4.5158 

-152.3782 

10.2705 

-2.8292 

-1.4805 

7.2623 

3 

1.4497 

1.1809 

5.1878 

2.8308 

-143.7239 

6.3551 

-1.7156 

-1.9421 

4.4938 

A 

1.4933 

1.2389 

4.0864 

2.0260 

-144.5321 

5.0624 

-1.1667 

-.8312 

3.5797 

5 

-.537V 

1.3117 

3.4377 

1.5325 

-140.0450 

4.3698 

-.8307 

-.6959 

3.0899 

6 

1*58*7 

1.3669 

3. C 636 

1.2386 

-133.713J 

3.9452 

-.6053 

-.6330 

2.7897 

7 

1.6243 

1.44C7 

2.8129 

1.0351 

-126.8436 

3.6841 

-.4389 

-.5857 

2.6051 

8 

1.666f 

1 .i 177 

2.6449 

.3961 

-119.9156 

3.5105 

-.3064 

-.5547 

2.4823 

9 

1.7117 

1 .cCtfc 

2.5229 

• 7961 

-113.1683 

3.3967 

-.1941 

-.5284 

2.4018 

10 

1.7553 

1.6646 

2.4355 

.7314 

-133.3643 

3.3254 

-.0933 

-.5087 

2.3514 

11 

1.799L 

1 .7661 

2.2740 

.6879 

-90.7821 

3.2776 

-.0066 

-.4863 

2.317b 

12 

1.0427 

-.8e49 

2.3316 

• 6889 

-79.5318 

3.2598 

• 0668 

-.4790 

2.3051 

13 

1.8653 

2.0C17 

2.3.90 

.5939 

-70.2615 

3.2610 

.1657 

-.4618 

2.3059 

14 

1.935C 

2.122C 

2.2988 

• 7336 

-61.0233 

3.2769 

.2513 

-.4538 

2.3172 

15 

1.9737 

2.2513 

2.3030 

• 7896 

-53.1616 

3*3101 

.3343 

-.4463 

2.3406 

16 

2.C173 

2.39C2 

2. 2157 

.9459 

-45.9916 

3.3557 

.4156 

-.4302 

2.3729 

17 

2.0610 

2.5394 

2.313d 

.9237 

-34.5445 

3.4136 

• 5036 

-.4158 

2.4138 

18 

2.1047 

2.6995 

2.3259 

1.0036 

-33.3933 

3.4824 

.5891 

-.3957 

2.4625 

19 

2 . 1 43 3 

2.6696 

2.3347 

1.3809 

-26.6000 

3.5617 

.6710 

-.3659 

2.5185 


ADDITIONAL 

. AXIAL 

LOCATION NO. 2, 

AT X/D » 

-.3634 






I 

R/D 

B/bINF 

8/BINF 

B/8INF 

b-angle 

B/BINF 

8X/BINF 

8T/8INF 

8Z/BINF 



(PARALLEL) 

(PERP) 

1 JN-PLANE) 

(IN-PLANE ) 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

1 

1.3751 

1 • f 526 








2 

1. 41-54 

1 .1100 

7.2744 

4* 3676 

-153**331 

9.8632 

-2.7526 

-1.4105 

6.9743 

3 

1.4636 

1.1635 

5.C6;& 

2.7608 

-149.3021 

6.1736 

-1.6786 

-.9966 

4.3654 

4 

1.5C79 

1.2222 

3.99*0 

1.9699 

-144.8720 

4.9481 

-1.1392 

-•0»15 

3.4988 

5 

1.5522 

1.2848 

3.3557 

1.4868 

-140.3229 

4.2715 

-.8092 

-.6713 

3.0204 

6 

1.5964 

1.3529 

2.9944 

1*2*21 

-133.6909 

3.8695 

-.5872 

-.6146 

2.7361 

7 

1.6437 

1.425t~ 

2.7529 


-120.6392 

3.6204 

-.4241 

-.5702 

2.5600 

8 

1.685L 

2 .5<. 35 

2.5917 

. 8 ?** 7 

-113.5724 

3.4564 

-.2945 

-.5407 

2.4440 

9 

i. 729 2 

1.5874 

2.4751 

• 7765 

-1 j9. 3692 

3.3485 

-.1621 

-.5180 

2.3678 

10 

1.7735 

1.6777 

2.3925 

• 7151 

-99.4893 

3.2837 

-.0834 

-.4987 

2.3219 

11 

1.817o 

1.7752 

2.335 U 

.5785 

-89.4659 

3.2405 

• 0045 

-.4798 

2.2913 

22 

2.862C 

1.6604 

2.2963 

.6 8Gb 

-78.3673 

3.2273 

.0971 

-.4714 

2.2821 

13 

1.9^63 

1.9939 

2*2769 

.6918 

-08.9646 

3.2320' 

.1754 

-.4566 

2.2834 

14 

1.9596 

2. lltu 

2.2694 

.7324 

-59.9162 

3.2526 

• 2596 

-.4401 

2.3000 

15 

1.9948 

2.2474 

2.2761 

• 7910 

-52.1266 

3.2894 

.3434 

-.4415 

2.3259 

16 

2.0392 

2.3687 

2.28t9 

.8514 

-45.C421 

3. 3390 

.4254 

-.4260 

2.3610 

17 

2. 0634 

2.54C5 

2.2911 

.9297 

-33.7107 

3.4CG4 

.5129 

-.4111 

2.4045 

18 

£.1276 

2 .7) 36 

2.3C-51 

1.0106 

-33.1342 

3.4733 

• 5984 

-.3906 

2.4560 

19 

2.1719 

2.e7e9 

2.3163 

1.3086 

-27.92B9 

3.5555 

.6801 

-.3605 

2.5141 

ADDITIONAL AXIAL 

LOCATION NO. 45 » 

AT X/0 • 

— 4 • 0009 






I 

R/0 

8/8INF 

3/8INF 

8/8 INF 

B-AN6LE 

8/BINF 

8X/BINF 

8Y/8INF 

SZ/8INF 



(PARALLEL) 

(PER/) 

(IN-PLANE) 

(IN-PLANE) 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

1 

1.7691 

.2912 








2 

1.917C 

• 4C08 

2.4912 

1.4725 

-178.8044 

2.6546 

-1.C41C 

-.3217 

1.8771 

3 

Z.W&49 

• 43C 4 

1.2732 

.5957 

-176.6662 

1.7587 

-.4206 

-.3230 

1.2436 

4 

2.2128 

• «»!>64 

.9545 

.3767 

-147.7656 

1.4419 

-.2253 

-.1421 

1.0196 

5 

c • 3 60 6 

.4662 

.0351 

. 3236 

-121.7165 

1.3097 

-.1203 

-.1947 

• 9261 

6 

2.5(35 

.5276 

• 7866 

.3235 

-103.3526 

1.2566 

-.0528 

-.2226 

• 8885 

7 

2.o564 

. 57c i 

.7719 

.3391 

-90.C363 

1.2483 

-.CC02 

-.2398 

.8827 

8 

2.0.-4 3 

.6371 

.7437 

.3869 

-79.6431 

1.2757 

.0492 

-.2691 

.9021 

9 

2.9522 

.7115 

.8152 

.4227 

-71.5340 

1.3309 

.0947 

-.2835 

.9411 

19 

S.lwit 

.6023 

.8664 

.4682 

-64.7550 

1.4128 

.1472 

-.3123 

• 9991* 

11 

3.2474 

.9 

.9*53 

• 539* 

-53.1646 

1.5189 

.2010 

-.3238 

1.0740 

12 

2.3958 

1.C432 

.9313 

.5919 

-50.8015 

1.6429 

• 2641 

-.3238 

1.1617 

13 

3.5437 

1.1987 

.9712 

.6677 

-44.&330 

1.7834 

.3394 

-.3281 

1.2610 

14 

3.6916 

1.38C 7 

1.C212 

.7664 

-33.0186 

1.9401 

.4270 

-.3338 

1.3718 

15 

3.8394 

1.590-7 

1.0736 

• 8769 

-32.367* 

2.1101 

.5237 

-.3319 

1.4921 

16 

3.9873 

1.8312 

1.1344 

1.0192 

-27.6554 

2.2924 

• 6383 

-.3345 

1.6210 

17 

4.1352 

2 .L 971 

l.loll 

. 1.1505 

-22.7621 

2.4700 

.7502 

-.3148 

1.7522 

18 

4.2831 

2.2962 

1.2699 

1.3244 

-19.5910 

2.6749 

.8823 

-.3140 

1.8914 

19 

4.431C 

2 . 7( 52 

1.3554 

1.5093 

-16.9048 

2.8557 

1.0211 

-.31*3 

2.0193 

ADDITIONAL AXIAL 

LOCATION NO. 4 6* 

AT X/0 ■ 

-4.9175 






I 

ft/ D 

8/&INF 

a/eiNF 

8/3INF 

8-ANGLE 

B/BINF 

BX/8INF 

ar/BiNF 

8Z/BXNF 



(Parallel) 

(PERP) 

(In-PLANE) 

(IN-PLANE) 

(NORMAL) 

(RESULTANT) 

(RESULTANT) 

(RESULTANT) 

1 

1.771C 

.3865 








2 

1.9237 

.4039 

2.5379 

1.5090 

179.8520 

2.5828 

-1.0670 

.CJ28 

1.8263 

3 

2.0765 

.4250 

1.2651 

.5936 

-178.9165 

1.7169 

-.4196 

-.0*79 

1.2140 

4 

2.2293 

• 45C7 

-.9682 

• 3883 

-149.3779 

1.4095 

-.2363 

-.1398 

.9967 

5 

2.362C 

• 4619 

.6576 

.3409 

-123.6337 

1.2818 

-.1333 

-.2008 

• 9064 

6 

2.5346 

. 5203 

• 6(6* 

.3377 

-105.4436 

1.2307 

-•C 636 

-.2302 

• 8702 

7 

2.6075 

.5662 

.7936 

.3545 

-92.3712 

1.2231 

-.0104 

-.2505 

• 6649 

8 

2.6433 

• c2?6 

• 8141 

• 4012 

-81.7914 

1. 2507 

•0405 

-.2808 

• 8844 

9 

2.9930 

• 7o:- 1 

.6345 

.4358 

-73.6228 

1.3066 

• C869 

-.2957 

.9239 

10 

3.145b 

.79i 9 

.8649 

.4993 

-66.7349 

1.3890 

.1394 

-.3243 

• 9822 

11 

3.2985 

.8948 

.9212 

.5478 

-59.9554 

1.4962 

.1939 

-.3353 

1.0580 

12 

3.4513 

1.1 3t 4 

.9461 

.5900 

-52.4654 

1.6196 

.2576 

-.3353 

1.1452 

13 

3.6040 

1.1651 

.9923 

. 6740 

-46.0303 

1.7603 

• 3309 

-.3430 

1.2447 

14 

3.7566 

1.3667 

1 .0 433 

.7721 

-39.8629 

1.9170 

.4191 

-.3500 

1.3555 

15 

3.9C95 

1.5761 

1.0859 

.6805 

-33.7606 

2.0868 

.5176 

-.3460 

1.4756 

16 

4.0623 

1 . site 

1.1502 

1.0193 

-28.9304 

2.2694 

• 6306 

-.3487 

1.6047 

17 

4.2151 

2.t.di 7 

1.1860 

1.1500 

-23.5757 

2.4542 

.7453 

-.3252 

1.7354 

18 

4.3678 

2.3765 

1.2617 

1.3237 

-20.5710 

2.6500 

.6763 

-.3269 

1.8738 

19 

4.5206 

2.6851 

1.3340 

1.5138 

-17.0815 

2.8289 

1.C211 

-.3138 

2.0003 





Figure A.2(j) - continued 
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APPENDIX A 


VcLDC I T7 AND TEMPERATURE CONTOURS 
12 VELOCITT ( Tfc NPER AT JRE ) CONTOUR LINES FOUND 

11 POINTS IN CONTOUR LINE Or tf/VINF - ,190# T/TINF ■ 22.120 


X/D R/0 


.9853 

. 19C 1 

1.TC19 

• L 6 32 

1.0109 

• 17e3 

1.C186 

.1741 

1 .*'35 3 

. io( 4 

1.0527 

• 14C3 

1.0707 

* *133 

1.J712 

• 1126 

1 *6 68 6 

.0614 

1.1933 

."363 

1.10 1C 

C.OCCC 


21 POINTS IN CONTOUR LINE OF V/ViNF • .250# T/TINF ■ 21.A60 


X/D 

.944 7 
.96 LA 
.9781 
.9961 
1.-1242 
1.0326 
1.9453 
1. >516 
1.0704 
l..>697 
1.11.95 
1.1151 
1.129 a 
1.1496 
2 .261 4 
1.1705 
1.1906 
l.:923 
1.2631 
1 • 2v3 6 


U/0 

.3697 

• 3te9 

• 3666 
. 3617 
.3552 

• 246 1 
.3396 
.3356 
.3213 

• 3k*4 
.2542 

• 278C 
.2577 
.2253 
.21 48 

a 1611 

• 1251 
.1221 
.0420 


OcNSITT CONTOURS 

lv DENS ITT CONTOUR LINES FOUND 


9 POINTS IN CiiNlJUfc LINE Dr RHOMHOJNF ■ 4.230 


X/D R/D 


• 993C 

.1282 

1.3C36 

• ;z;; 

1.0247 

.1171 

l.L-34 6 

• lv£6 

1.0417 

..•925 

1.0597 

.-376 

1.0598 

♦ 0*7i 

1.0599 

.0242 

1.3599 

•; ■ -'c C-C 


Figure A.2(k) - continued 
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AFPSHDIX A 


tr* POINTS IN CONTOUR LINE OF RHO/RHOINF - ♦•000 


tf 0 

.9983 
•9ei7 
.9759 
• 9 895 
1.0C39 
1.C189 
1.3399 
1.0508 
1 .36*1 
1.0667 
1.1071 
1 *1161 
1.1298 

1.1558 

1.1559 
1.1795 
1 • 1 87 ♦ 
1.1939 
1 .?tC 1 
1 .2005 


P/0 

.3572 

• 3o73 
.3776 
. 3661 

. 3933 

• 3966 

• ♦*.16 
.♦CI6 

.3977 

• 386 6 
.3722 
.3626 

• 3906 
.2882 

• 2eei 

.2i a: 

• U91 
.1255 
.0919 

C . JC cc 


MAGNETIC FIELD CONTOURS 


11 MAGNETIC FIELD CONTOUR LINES FOUND 

(FOR COMPONENT ALONG FIELD LINES PARALLEL TO FLOW IN FR6ESTAEAN) 


13 POINTS IN CONTOUR LIN: OF 8/8INF (PARALLEL) • .500 


X/D 

.9789 
.9939 
1.0323 
1.0296 
1.0973 
2.75>9 
1.0699 
1.0831 
:.>9&6 
l.lf 16 
1.1176 
1.1202 
1.12ZC 


R/D 

.2283 

• 2726 

• 23 99 
.2833 

• 1692 
.1859 
.1689 
.1905 
.1253 

• 39e3 

• i:9t 

.Jlt8 

3.0CCC 


19 POINTS IN CONTOUR LlNt JF 8/8INF (PARALLEL) - .750 


X/D 

.9513 
.9669 
• 9 86C 

1.0193 
2.0122 
1.0222 
1.0903 
1.C589 
1.3 78 C 
1.0911 
3.3973 
1.1166 
1.1369 
1.1927 
3.156* 
1.1799 
1.1767 
1.1 683 
1.1892 


A/0 

• 3969 
.3 937 
.3391 
.3326 

• 3289 
.3292 
.3236 
.3CC9 
.2895 
.2721 
.2651 

• 2 *C 6 

.2116 

• 2115 
.1691 

• 1239 
.1108 
.0*25 

G.OCCC 


Figure A. 2(1) - continued 
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APPENDIX A 


12 MAGNETIC F1ELL CONTOUR LINES FOJNO 

(FOR COMPONENT ALONG FIELO LINES PEAPENOICULAR TO FLQ» IN FREESTREAM) 

NS POINTS IN CONTOUR LINE OF 3/EINr (PERPENOICULaRI - R.Jili 


X/D 

R/D 

-1.8975 

1.7899 

-1.0527 

1.7957 

-1.7055 

1.8342 

-1 .6156 

2.9226 

-1.5257 

l.o283 

-1.5C27 

1.3293 

—1.4359 

1.3246 

-1.3459 

1.3394 

-1.2560 

1.8414 

-1.166? 

i.84f 7 

-1.1279 

1.8396 

-1.0762 

1.8425 

-.9863 

6.8 446 

-•34o4 

1.0430 

-.6542 

1.84C9 

-.81 65 

2.3452 

-•753C 

l.o466 

-.6996 

1.3499 

-.6814 

1.65C' 

—.6461 

2.0559 

-•5926 

2.062: 

-.5626 

1.8652 

-.5392 

1.0726 « 

-.4657 

1.8677 

-.4638 

1.3882 

-.4322 

I.924& 

-.4277 

2.9265 

-.3885 

1.9497 

-.3786 

1.9667 

-.3651 

1.9690 

-.3555 

2.0345 

-.3fc2e 

2.v697 

-.3788 

2.1363 

-.3845 

2.1524 

-.4257 

2.2281 ; 

-.4322 

2.2404 

-.4 70 6 

2.3174 

-.4857 

2.3311 

-.S3C7 

2.3973 

-.5392 

2.4C9J 

-.5926 

2.4803 

- .61.96 

2.51.17 

-.6461 

2.5431 

-.6996 

2.61 K 

-.7337 

2.6367 

POINTS IN CONTOUR LlNt 

OF d/BlNF (Pc0P£NDICULAR) ■ 2.000 


X/D 

R/0 

-3.6973 

1 .0865 
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Figure A. 3. - Plot output for sample case 
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Table 1 (Concluded) 
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Figure 1 . - Comparison of the equatorial and principal meridian 
traces of the magnetosphere boundary as provided by the 
simplified theory of equation 16 . 
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Figure 5. - Comparison of flow properties predicted by the present implicit method 
with other techniques and experiment for supersonic flow past a sphere; 
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(a) Shock shape and sonic line location. 



x/D 

(b) Density distribution. 

Figure 6. - Comparison of implicit and inverse methods for 
shock shape and sonic line location, and density dis- 
tribution along bow shock and magnetosphere boundary 
for =8,7= 5/3 flow past the rotated 
equatorial trace of the magnetopause. 
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Figure 7 . - Bow wave and sor 
for various supersonic f 
rotated equatorial tr 
magnetopause with K 


3.0 
















(a) (Magnetopause equatorial trace 


(b) Ionopause trae<e 


(c) Ionopause trace - H/R 0 ' = .3 - ~ 

Figure 11. - Variation of shock stand-off distance 
with oncoming Mach number andtrratio of specific 
. heats for various magneto/ionopause traces as 
determined by the present implicit procedure. 
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Figure 12. * shock shapes for various, supersonic flows past 
the rotated equatorial. trace of the magnetopause; 
combined near, (blunt body) and far 
(marching) solutions. 
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Figure 13. - Shock shape for = 8, Y = 5/3 flow past an ionopause shape with 
H/R_ = 0.1; combined near (blunt body) and far (marching) field solutions. 












Figure 14 . - Continued 
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(a) Perpendicular component ( B/B „,) A 


Figure 15. - Contours and field line locations of the in-plane magnetic 
field components (B/B^Jx and (B/B OT ) ,, for M,,, = 8 and , Y = 5/3 flow 
past the rotated equatorial trace of the magnetopause. 
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